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Ce cours de resistance des materiaux a pour objectif d'approfondir la mecanique des 
solides elastiques, puis a partir de la mecanique des milieux continus, nous introduirons la 
theorie des poutres. Dans une premiere partie, nous etudierons la demarche qui nous permet 
l'etablissement des equations de la theorie des poutres (une demarche similaire pourrait etre 
utilisee pour les plaques et coques). Dans une deuxieme partie, nous nous attacherons a 
exposer les outils classiques de la theorie des poutres: etude de cas simples, methodes 
energetiques , . . . etc . . . 

La resistance des materiaux est un outil indispensable a toute modelisation en calcul des 
structures. Meme si d'autres methodes (par exemple les elements finis) sont en general 
utilisees, un calcul rapide de RDM permet de verifier les ordres de grandeur et de juger de 
l'opportunite d'utiliser d'autres methodes plus complexes. 

Ce polycopie est en perpetuel correction (quand j'en prends le temps). C'est pourquoi, je 
serai reconnaissant aux etudiants de m'exposer toute suggestion susceptible d'en ameliorer le 
contenu. 
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RAPPELS DE MECANIQUE DES 
MILIEUX CONTINUS 



I - Cinematique 



I - 1 Configuration, mouvement, deplacement, ... 

L'espace physique est rapporte a un repere orthonorme direct (O, ei, Si, e?)- L'ensemble 
des particules ou points materiels constituant le milieu continu etudie, occupe a chaque instant 
t, un ensemble de positions dans l'espace: c'est la configuration du systeme a l'instant t, note 
Q(t) (d'interieur Q(t) et de frontiere 3Q(t)). 

On introduit aussi la notion de configuration de reference: c'est la configuration 
particuliere du systeme a un instant to fixe. Souvent on prendraQ = Q(0), et on parlera alors 

de configuration initiale. 

Toute particule Mq de Q, est reperee par son vecteur position X(t) dans la configuration 

de reference. Toute particule M de Q.(t) est reperee par son vecteur position x(t) dans la 
configuration actuelle (a l'instant t). 

n n _ 

~ x(X,t) 



o 



ei 




Configuration cCe reference 
a ['instant to 



Configuration actuetk 
a ['instant t 



Figure 1 

La position de chaque particule M sera done determinee si on connait sa position dans la 
configuration de reference et une fonction $ telle que: 

x(t) = cl>(x,t) (1) 

4> definit le mouvement par rapport a (0,ei,e*2, e?) pjj re q Ue \ e m iij eu es t continu, c'est 
dire que <J> est une fonction continue et biunivoque de X. 

X et t definissent les variables de Lagrange 

x et t definissent les variables d'Euler 
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- 8 - Rappel de MMC 

Le emplacement par rapport a la configuration Q. , a l'instant t, de la particule Mo est le 
vecteur 

u(X,t) = x(X,t)-X (2) 



I - 2 Deformation 



Considerons deux particules voisines X et X+dX. A l'instant t ces particules occupent la 
position x et x+dx avec dx(t) = 4>(x + dX, t)- <I>(x, t) 



Par definition du gradient on ecrit 



<J>(x + dX, t) = 4>(x, t)+ —fx, tjdX + 

oX 



dX 



2\ 



Soit 



dx = F(x,t)dX avec F(x,t)= ^-(x,t) 

oX 



(3) 



F est une application lineaire qui fait passer de l'espace vectoriel dans lequel peut varier 
dX dans l'espace vectoriel ou varie a priori dx . Cette application lineaire, appelee tenseur 
gradient, permet done le passage de la configuration Q a la configuration Q(t) . 

En notation indicielle, 

dx 1 



Fij "ax, 



3xi 



ax ; 



F = 



3x 1 



3x 1 



axj 


ax 2 


ax 3 




dx 2 


dx 2 


dx 1 


ax 2 


ax 3 


a x 3 


a x 3 


a x 3 


axj 


ax 2 


ax 3 



Remarques: 

* Transformation d'un element de volume dV dans Q. en un element de volume dv dans Q(t) . 

dv = det(F}lV 

* Transformation d'un element de surface N dS dans Q. en un element de surface dans Q(t) . 

iids = detFF T NdS 



Le tenseur gradient decrit la transformation locale au voisinage d'une particule donnee. 
Afin de rendre compte des deformations, e'est a dire des changements de forme autour de 
cette particule, on s'interesse a revolution du produit scalaire de deux vecteurs materiels pris 
respectivement dans les deux configurations Q. et Ci(t). 

Considerons trois particules voisines X, X+dX, X+dX'. Apres deformations, elles occupent 
dans Q(t) les positions respectives x, x+dx, x+dx'. 
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dl ■ dx = (f(X, t) dx)- (f(X', t) dX')= 



3xk 

v ax t 



dXj 



34 

, ax- 

V J 



dX- 



d'ou sa variation autour de la transformation 



dx-dx'-dX-dX' = 



. axi ax ■ lj 



dXidXj 



soit 



dxdx'-dXdX' = 2dX£dX' 
en po sant 




L' application lineaire £ est appelee tenseur des deformations. Cette application est 
symetrique mais depend bien sur de la base (O, ei, 82, 83) initialement choisie. 



Autre, ecriture: 
D'apres (2) et (3) 



F(X,t) = ||(X,t) = l + J|(X,t) 



soit 



£ = 



au . au /xr . au . au . 

-<X,t) + - (X,t) + - (X,.)-(X,,) 



(5) 



ou encore en notation indicielle 



1 



^au L+ au i+ au k au t ^ 
axj axi axj axj 



1-3 Cas des petites perturbations 



Cette hypothese correspond au cas ou ||u(X, t)|| et 



sont petits. 



En re prenant (5) et en ne retenant que les termes d'ordre 1, on obtient: 



1 

£hpp — — 
2 



^-(X,t) + ^ (X,t) 

dx dx 



(6) 



ou encore en notation indicielle 



1 

e ii =- 
1J 2 



auj 3uj 
axj axj 
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I - 4 Conditions de compatibility 



A tout emplacement u on fait correspondre une deformation £ . On peut aussi se poser le 
probleme inverse. Ce probleme est dit 'probleme de compatibilite geometrique d'un champ de 
deformation', ou encore 'probleme d'integrabilite d'un champ de deformation' . 

Les conditions de compatibilite peuvent etre etablies dans le cas general, cependant nous 
ne les etablirons que dans le cas des petites perturbations. 

Decomposons maintenant le gradient des deplacements en une partie symetrique £ et une 
partie antisymetrique CO . 



||(X,t) = £(X,t) + 00(X,t) 
du . 3u . 



©ii =- 
1J 2 



3uj 3"j 

axj ax t 



On a 



GO 



ij,k 



^ki,j ^jk,i 



soit en derivant une nouvelle fois 



ij.ki = ro ij,ik v dans { i' 2 * 3 ) 



Vi,j,k,l£y kl +£uij — e ik,ji — e ji,ik -0 



(7) 



Reciproquement, si £ verifie (7), alors les formes differentielles 
d0) ij =( e ki,j-£jk,i) dx k 

sont exactes; elles permettent done de construire le champ CO de tenseur antisymetrique. 
On verifie ensuite que les formes differentielles 
du ; =((0 ik +£ ik )dx k 

sont exactes, d'ou la possibilite de construire un champ de deplacement u(X,t) defini dans Q. Q . 
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II - Sthenique 

II -1 Forces 

Elles resument les effets mecaniques, autres que cinematiques, exerces sur le milieu 
continu considere par le reste du domaine physique. Leur schematisation a chaque instant 
repose sur la definition d'un champ de vecteur 3>(x,t) et d'une mesure positive (0, definis sur 
la configuration actuelle Q(t) .<£(x,t) est une densite de force pour la mesure tt>. 

* Si (0 est une mesure de volume, alors <3>(x,t) est une force volumique (densite 
volumique de force) definie dans Q(t) de la configuration actuelle, par la fonction 

f: xen(t)->f(x,t)e 9\ 3 

* Si (0 est une mesure de surface, alors $(x,t) est une force surfacique (densite surfacique 
de force) definie sur 3Q F (t) de la configuration actuelle, par la fonction 

F: xe3Q F (t)^F(x,t)e 9x 3 

* ... etc ... 
Remarques: 

* Les forces sont definies sur la configuration actuelle. 

* A un instant donne et en un point donne x de d£l(t), on ne peut imposer a la fois le 
deplacement et la force !. Mais l'un des deux doit etre impose. On note 3Q F (t) la frontiere ou 
la force est imposee, et 3^u(t) la frontiere ou le deplacement est impose. Dans le cas des 
appuis mobiles, les composantes non imposees cinematiquement le sont pour les forces 

* Le monde exterieur au milieu considere doit, pour imposer le deplacement U(t) au bord 
d£2u(t), exercer des forces que nous noterons R(x,t) . Comme elles sont a priori inconnues, 
nous les appellerons reactions pour eviter de les confondre avec les autres forces qui, elles, 
sont donnees. 

II - 2 Contraintes 

II - 2.1 Notion de Vecteur-contrainte et tenseur des contraintes 

Soit un corps (C) en equilibre par application d'un 
sy steme d'actions mecaniques exterieures. Imagi-nons 
qu'une surface £ divise (C) en deux parties (1) et (2). La 
| partie (1) est en equilibre sous les actions mecaniques 
exterieures qui lui sont appliquees et les actions 
mecaniques exercees par la partie (2). Nous admettrons 
que sur chaque element de surface dE deE, (2) exerce sur 
(1) une force dF(x,t,h) 1/2 de densite superficielle 
T(x,t,n). 

dF(x,t,n)i/ 2 =T(x,t,n)d£ (8) 
T(x,t,h) est le vecteur contrainte au point x, relativement a la facette dE definie par son 
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vecteur normal n . 

La densite surfacique de forces exercees en x depend de x, t et aussi de l'orientation de la 
surface E au voisinage de x. Elle est lineairement dependante de h. On introduit alors 
l'application o telle que: 

T(x,t,n) = a(x,t)n (9) 

L'application o(x,t) s'appelle le tenseur des contraintes de Cauchy en x a l'instant t; il 
caracterise, dans la configuration actuelle, les efforts interieurs de cohesion exerces sur une 
partie du solide a travers l'element de surface ndZ. 

II - 2.2 Autre ecriture du tenseur des contraintes 

En utilisant la remarque du §1-2 pour exprimer ndZ en fonction de N dS, (8) devient: 
dF(x(X, t), t, n(N, t)) = n N(X)dS 

ou II est le tenseur 

fl(X,t): Ne9t 3 ^n(X,t,N) = fl(X,t)N e <K 3 
defini par 

n(X,t) = (detF)aF T 

Cette application lineaire n(X,t), definie pour X e Q, , s'appelle le premier tenseur des 
contraintes de Piola-Kirchoff en X a l'instant t; la composante IL; est la i ltaie composante du 
vecteur contrainte exercee sur la deformee d'une surface unite, normale a ej, de la 

configuration de reference. On prendra garde aufait que le tenseur II n'est pas symetrique. 

Si maintenant on cherche le vecteur "force de cohesion" dans la configuration de reference 

dF (x,t,N)=F 1 (X,t)dF(x(X,t),t,h(N,t))=SN(X)dS 
ou 3 est le tenseur defini par 

S = F 'n 

Cette application lineaire S(X,t), definie pour X eH , s'appelle le second tenseur des 

contraintes de Piola-Kirchoff en X a l'instant t. Attention, sa composante Sy n'est pas la 

jieme composante du vecteur contrainte exercee sur la deformee d'une surface unite, normale a 
e j , de la configuration de reference, mais seulement la i lhme composante de son transporte 

dans la configuration de reference. 

Selon le jeu d'ecriture adopte, on a done trois descriptions des contraintes: 



o = (d 



=\-l=T / =\-l 



^detFj IIF =|detFj FSF (10) 
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II -3 Equilibre 



II - 3.1 Le Principe des Puissances Virtuelles 

Pour schematiser les efforts mis en jeu, il est commode d'imaginer des mouvements fictifs 
(ou virtuels) et d'analyser le travail ou la puissance qui en resulte. Par exemple, pour evaluer 
les forces de gravite agissant sur un objet, on peut imaginer de le soulever (mouvement virtuel 
de bas en haut). 

II-3.1.1 Le Principe des Puissances Virtuelles (Germain 1972) 
Un milieu materiel etant isole, on peut distinguer les actions exterieures qui agissent sur le 
milieu, des actions interieures qui representent les liaisons existant entre toutes les parties du 
milieu. 

Axiome d'objectivite 

La puissance virtuelle des efforts interieurs associee a tout mouvement rigidifiant est nulle. 

Axiome d' equilibre 

Pour tout milieu materiel repere dans un referentiel absolu, a chaque instant et pour tout 
mouvement virtuel, la puissance virtuelle des quantites d 'acceleration Y\ a est egale a la 
somme des puissances virtuelles des efforts interieurs Yli et des efforts exterieurs Y\ e . 



II-3.1.3 Puissance virtuelle des efforts interieurs 
Pour determiner la puissance virtuelle des efforts interieurs nous ferons les hypotheses 
suivantes: 

* IJj admet une densite volumique pf 



* TIi est en chaque point une forme lineaire des valeurs en ce point de Sv et de ses 
derivees premieres: 



II-3.1.2 Position du probleme 




Soit un milieu continu Q.(t) d'interieur Q.(t) et de 
frontiere dQ.(t). Isolons maintenant un domaine Z(t) 
de frontiere 3E(t) interieur a £2(t), et soit n la 
normale en un point de dZ (t). A un instant t fixe, 
un mouvement virtuel defini par une vitesse 
virtuelle Sv est applique a E(t). Cette vitesse est 
supposee continue et continument derivable sur 
E(t). 



n i =!!!zPi dx 
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En decomposant le gradient des vitesses virtuelles en une partie symetrique 8D et une 
partie antisymetrique 8W , 



sw=- 



8D=- 



asv asv 

dx dx 



asv asv 
— + — 

dx dx 



8W, 



1 



1J 2 
1J 2 



la densite volumique des efforts interieurs devient: 
Pi = Aj 8vi + By 8wij - Oy 8Dij 



aSvj 35vj 
ax, dx { 

asv: asv^ 

L + — 

axj ax ; 
(ii) 



Le premier axiome du principe des puissances virtuelles impose que pour tout mouvement 
de solide rigide la puissance des efforts interieurs soit nulle. D'ou: 

- Soit un mouvement de translation: 8v ^ 0, 8W=0 et 8D=0 
alors 

n i =MIi Pi dx= !!J£ A ■ 8vdx=0 V X dans Q 

soit A ■ 8v = V 8v, ou encore A = 

- Soit un mouvement de rotation: 8v = 0, 8W^0 et 8D=0 

alors 

nj = Pi dx = B:8W dx = V £ dans Q 

soit B : 8W=0V8W , ou encore 1=0. 



Done en definitive: 

I!; =-||| £ o:8Ddx (12) 



On peut montrer que le tenseur o introduit ici correspond bien au tenseur des contraintes 
de Cauchy exprime au §11-2.1. 



II-3.1.4 Puissance virtuelle des efforts exterieurs 
Les efforts exterieurs comprennent 

- des efforts exerces a distance par des systemes exterieurs a £1, supposes definis par une 
densite volumique de forces f , 

- des efforts de cohesion schematises par une densite surfacique de force T sur 

n e =fffcf •Svdx+Jj 3E f -Svdx (13) 
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II-3.1.5 Puissance virtuelle des quantites d'acceleration 
Si y est l'acceleration et p la masse volumique de chacun des points de alors 



n a =|jjvPY-5vdx (14) 



II - 3.2 Application du Principe des Puissances Virtuelles 

En application du Principe des Puissances Virtuelles on obtient: 

- JJk o:5Ddx+|j| E f ■ 8vdx+|| as f ■ 6vdx=|jj £ py ■ 8vdx (15) 

Pour exploiter le fait que (15) est verifie pour tout mouvement virtuel, nous allons faire 
apparaitre 8v dans chacun des termes. 

En appliquant le theoreme de la divergence, le premier terme devient: 

= = = 3§v = = 

- a:8Ddx=- ||} E a:— — dx=-j] a£ a • 8v • ndx+[fj* E div x o • 8vdx 

OX 

Soit: 

U aE (f-a ■ n) 8vdx+|jj E (f +div x o-py)- 8vdx = 0V8v 
Ou encore 



f +div _ o=pydansE 

, = (16) 
T=o ■ n sur3Z 



II - 3.3 Equilibre 

En considerant les developpements du paragraphe precedent et en se ramenant au domaine 
Q(t), nous pouvons done ecrire les equations d'equilibre d'un solide soumis a un champ de 
forces exterieures f dans O(t), a un champ de forces exterieures F e sur 9Q F (t) et a un 
deplacement impose U e sur 3Qu(t)- 

Dans la configuration actuelle: 

f (x, t)+div x a(x, t)=0VxeQ(t) 

= f , f , {F e (x,t)Vx G 3^ F (t) 
o(x,t)-n(x,tH - 

[R(x,t)Vxe3nu(t) 

Dans la configuration de reference: 
De meme, si on note fg, R et F les densites volumiques et surfaciques de forces 

mesurees dans la configuration de reference: 

f (X,t)+div x rT(X,t)=OVxG^ (19) 



(17) 
(18) 
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|F (x,t)Vxex" 1 (3a F (t),t) 

n(x,t)-N(x,tH 

R (x,t)Vxe3a 0U 



Rappel de MMC 
(20) 



II - 3.4 Cas des petites perturbations 

Reprenons (17), en l'exprimant en fonction de X 

fi(x(x,t),t)+^i(x(x,t),t)=ovx(x,t)Gn(t) 

d Xj 

f i (x(x, t), t)f|^- (x, t)f^ (x, t)=o vxg a 

dX k dXj 

Or x(X,t) = X+ u(X,t) soit -^-(X, t)=l+ -^-(X, t) 

oX oX 

On peut done ecrire l'equation d'equilibre sous la forme 



fi(x(X,t),t)+!^(X,t) 
dX k 



7 ^ X 

1+ S (X -° 



=0VXeQ 



kj 



Sous l'hypothese des petites perturbations, on peut alors ecrire: 



1+ S (X -° 



soit 



30;; I 
fi(x(X,t),t) + — i (X,t) 
dX k 



5 ^k (X ,t) 

Jk ax. 



1 



= VXeQ f 



Enfin, en ne retenant que les termes d'ordre 0, et apres avoir effectue un developpement de 
fi au voisinage de X, on obtient: 

fi(X,t)f^-(X,t)=0VXen 



soit 



f (X, t)+div x o(X, t)=0VxeQ 



(21) 



Le raisonnement qui a permis de remplacer f(x(X,t),t) par f(X,t), permet aussi de 
remplacer F e (x(X,t),t) par F e (X,t) et R(x(X,t),t) par R(X,t). Done, comme condition sur la 

frontiere on obtient: 



fF e (X,t)VXGdQ 0F 
a(X,t)-N(X,tH 


(22) 


[R(X, t)VXe3Q 0U 





I.S.I.T.V. 



Resistance des Materiaux 



-17- 



II - 3.5 Autre approche 

On utilise la loi fondamentale de la dynamique qui stipule que : 

le torseur dynamique, qui est la derivee temporelle du torseur cinematique est egal au 
torseur des actions exterieures. 

Ce qui se met en equations sous la forme suivante, en l'appliquant au domaine d 'etude E : 



-JJJ E pvdx = JJ ffi fdx + JJJ E fdx 
dt 



(23) 



— JJLOM Apv dx = JLOM aT dx + JJLOM a f dx 
dt 



Par application des proprietes de la derivee particulate on peut ecrire pour tout vecteur b : 



dt 



3pb 

at 



+ div(pb®v) 



dx 



puis en utilisant le fait que pour tout vecteurs a et b 
div(a ® b) = Va • b + a divb 

on peut developper sous la forme 



dt 

d_ 
dt 



JfJ E pbdx= E 



jjj zP bdx= in, 



apb 
at 



+ Vpb ■ v + pb divv 



dx 



J 



p ^ + b ^ + p Vb ■ v + Vp ® b ■ v + pb divv 

at at 



A|}} E pbdx= nt 

dt 



v 
f / 

P 

V v 



dx 



5b V7U - 

— + Vbv 

at 



+ b 



J 



— + Vp ■ v + p divv 

at 



dx 



Par definition de la derivee particulaire 



dt 



f z pbdx= JJt 



db rfdp 



p — + b 
dt 



at 



+ divpv 



dx 



(24) 



Done, grace a la conservation de la masse 
— + divpv = — + p divv = 

at dt 

et la definition du vecteur contrainte (8), la premiere equation de (23) devient : 

^ JJtpv dx = {{{,p ^ dx = f ,pY dx = jj 3z 5n dx + JJJ J dx 
dt dt 

et par application du theoreme de la divergence 
JIJzPY dx = jjj z diva dx + \\\J dx 

On r etrouve done bien la forme locale de la conservation de la quantite de mouvement : 



divo + f =py 



(25) 
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II - 4 Quelques proprietes du tenseur des contraintes 



Le tenseur des contraintes est un tenseur symetrique. Dans tous les developpements a 
venir, nous nous placerons dans le cas des petites perturbations pour un solide en equilibre. 
En consequence, nous omettrons les variables x et t. 

II - 4.1 Contrainte normale et contrainte tangentielle 



Considerons une facette de normale n . Tout naturellement, le 
,■> vecteur contrainte T(n) peut etre decompose en une 
composante normale G n et une composante tangentielle T. 

o n =T(5)-5=n-o-5 
et 




11x11=^(0-5)* -(s-o-s) 
On dira que a n est positive en traction et negative en compression. 



II - 4.2 Directions principales, contraintes principales 

La matrice representant le tenseur des contraintes est symetrique, elle est done 
diagonalisable. Les valeurs propres sont reelles et appelees contraintes principales (Oi, On, 
Oni)- Les vecteurs propres, orthogonaux deux a deux, sont les directions principales 
(5 l5 n n , n m ). On a done: 

o I= f (5 T ) ■ n I ,o n =f (n n ) ■ 5 n ,o m =f (5 m ) ■ n m 



II - 4.3 Invariants 

Le tenseur des contraintes possede trois invariants definis mathematiquement comme les 
coefficients de l'equation caracteristique det^ - a l). C'est a dire les quantite scalaires: 



Tr(o) 
1 



Tr(o) z -Tr(o ) 



E m =Det(o) 

Exprimes en fonction des contraintes principales, on obtient 
Ei = Oi + o n + o m 
E n = Oi o n + o n o ra + o in &i 
Z m = Oi o n o m 



(27) 
(28) 

(29) 



II - 4.4 Cercles de Mohr 

Connaissant le tenseur des contraintes o, on se propose de determiner le domaine 
engendre par l'extremite du vecteur contrainte quand n varie. Par commodite, nous nous 



I.S.I.T.V. 



Resistance des Materiaux 



-19- 



pla£ons dans une base orthonormee dirigee suivant les directions principales de a . Soit 





V 









0" 




n=< 


n 2 


>eta= 





°2 





2 2 2 

avecn 1 +n 2 +n 3 =l 




. n 3, 










°3_ 





et 



T= 



n 2 °2 
n 3 ^3 



D'apres (23) 



o n = o T n, + o TT n 9 +G T 



•n - °I n l TU II U 2 _r ' J in n 3 
et d'apres (24) 

2 2 2 2 

x +o n = Oj n 1 + c 

Dans l'hypothese ou les contraintes principales sont distinctes, on obtient alors apres 



2 , ~ 2 2 
2 + °III n 3 



resolution du systeme: 



2 _x +(o n -an)(o n -o m ) 



n, =■ 



(o I -o II )(o I -o III ) 



n 2_ x +(o n -Oi)(o n -ani) 

n 2 — 

(^II-^l)^!- !!!) 

^2_ x 2 +(a n -Oi)(o n -o n ) 
(o m -Oi)(Oni-On) 

Si on ordonne les contraintes principales de telle sorte que Oi > On ^ Oni , alors 

x 2 +(o n -o„)(a n -o ni )>0 

x 2 +(o n -o I )(a n -a III )<0 

x 2 +(o n -o I )(o n -o II )>0 
ou encore 



x 2 + 



Ott+O 



II™ III 



x 2 + 



x 2 + 



2 / 



, n -u m 



°i+°m 



o T +o n 



< 



a T -o 



in 



o T -o n 



\ 2 



(30) 
(31) 
(32) 



j 



Dans le plan de Mohr, l'extremite du vecteur contrainte, d'apres (31), est done interieure au 
cercle centre sur 0o n d'abscisse (Oi + Oni)/2 et de rayon (Oi - Oni)/2. Par contre, d'apres (30) 
(res. (32)), l'extremite du vecteur contrainte est exterieure au cercle centre sur 0o n d'abscisses 
(On + <5in)/2 (resp.(Oi + o n )/2) et de rayon (o n - o m )/2 (resp.(Oi + a n )/2). 
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Description des Cercles principaux: 

Nous allons etudier la description du grand Cercle de Mohr. Les facettes concernees sont 
paralleles a la direction associee a la contrainte principale On- 

On constitue avec les directions I,III,II un triedre direct 
(Ojejjinjjin), la normale n de la facette evoluant dans le plan I 
III. 

Et on definit Tangle = (Ln), et le vecteur t tel que (n,t,II) soit 
direct. 

On a alors 

ii^osGej+SinGejjj 

et 

^OjCosGej+OjjjSinGejjj 

En utilisant les formules de changement de base de (0,ej,e in ,e n ) a (n , t ,11), on a done 
o I+ o m | a I -a m cQs2e 
2 2 

x= _^in Sin2 
2 

Lorsque la facette tourne autour de la direction de 
la contrainte principale Oji d'un angle donne, 

I'extremite du vecteur-contrainte tourne sur le cercle 
de Mohr d'un angle double dans le sens oppose 
(autour du centre du cercle). 
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III - Loi de Comportement pour les solides elastiques 

Pour determiner revolution d'un systeme deformable, nous avons deja determine les 
equations de la cinematique et de la sthenique. A ces equations, il est maintenant necessaire 
d'adjoindre une relation supplementaire reliant les efforts internes et les grandeurs 
cinematiques. Cette relation, appelee Loi de Comportement, depend du materiau considere. 
La construction d'une loi de comportement est basee sur des observations experimentales. 

Dans ce chapitre nous exposerons le modele de comportement des materiaux elastiques, 
sous l'hypothese des petites perturbations. 



Ill - 1 Approche experimentale: essai de traction 



So 



Lo 





Pour effectuer un essai de traction simple sur un 
metal, on utilise une eprouvette cylindrique 
caracterisee par: 

- des extremites surdimensionnees 

- des conges de raccordement (pour eviter les 
concentrations de contrainte) 

- une partie mediane cylindrique dans laquelle le 
champ de contrainte est suppose homogene, de 
traction simple parallelement a l'axe de l'eprouvette. 



L'essai de traction consiste a enregistrer revolution 
de l'allongement relatif de la longueur initiale Lo en 
fonction de la force de traction F, ou du rapport 
F/So, ou So represente l'aire initiale de la section de 
l'eprouvette. 

La figure ci-contre represente un tel enregistrement pour un 
acier inox. On remarque alors les proprietes suivantes: 

- Lediagramme est independantde la vitesse dechargement 

- La partie OA du diagramme est reversible. Si on charge 
jusqu'a un niveau inferieur a Go, alors la decharge decrit la 

meme courbe OA. 

- La partie reversible est lineaire 

- Si on effectue un chargement au dela du seuil Oo, puis 
une decharge, l'eprouvette presente une deformation 



permanente. 



I.S.I.T.V. 



-22- 



Rappel de MMC 



La partie reversible du diagramme de traction est, par definition, representative du 
comportement elastique du materiau. o est la limite initiale d'elasticite du materiau. La 

linearite du segment OA caracterise le comportement elastique lineaire du materiau. 



Ill - 2 Loi de comportement elastique lineaire (en HPP) 

III - 2.1 Forme generate 

A partir des observations experimentales on peut ecrire que les contraintes dependent 
lineaire ment des deformations. En l'absenced'effetsthermique et de contraintes initiales on a: 

o(x,t)=C(x):e(x,t) (33) 

C est un tenseur du quatrieme ordre, dont les composantes sont les coefficients d'elasticite 

du materiau. 

o ij (x,t)=C ijkl E kl (x,t) 

En utilisant les proprietes des tenseurs de contrainte et de deformation, on peut montrer que: 

Qjkl =_Cjikl Qjkl = Qjik Cijkl = Ckiij 

Le tenseur C, dont la matrice representative comporte 81 composantes, ne depend done 
plus que de 21 parametres independants. 



Ill - 2.2 Materiau elastique homogene isotrope 

Toutes les directions sont equivalentes, de telle sorte que la loi de comportement est 
invariante dans toute rotation de la configuration de reference. Ce modele s'applique a la 
plupart des materiaux: acier, beton, ... 

Si la configuration est libre de contraintes, alors la loi de comportement s'ecrit: 



o=XTr(8)l+2(iE (34) 



ou encore en notation indicielle 
Oij=^kkS ij +2u£ ij 

Les coefficients materiel X et (I, qui dependent de la particule consideree, sont appeles les 

coefficients de Lame. Leur expression en fonction du module d'Young E et du coefficient 

de Poisson v, est 

E n vE 
u= etA= ou 

P 2(l+v) (l+v)(l-2v) 

E = — ^etv = 

X + (I 2(X + |i) 

avec, en inversant (34) 
v = = l+v = 

8=—Tr(o)l+— o (36) 
E E 
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III - 2.3 Materiau elastique homogene orthotrope 

Le materiau possede trois directions privilegiees deux a deux orthogonales. La loi de 
comportement est invariante par les symetries par rapport aux plans orthogonaux construits a 
partir de ces directions. Dans ces materiaux, on peut classer les toles laminees, les composites 
tisses, le bois, certains betons structures, ... 

Dans ce cas on montre que la matrice de comportement est definie par 9 parametres 
independants. Dans le repere principal d'orthotropie, la loi se met sous la forme: 





1 


-V12 


-Vi 3 















eT 


Ei 


Ei 










Ell 




-v 2 i 
E 2 


1 

eT 


-v 23 
E 2 















£ 22 




-v 3 i 


-v 32 


1 















e 33 




E 3 


E 3 


eT 


1 






< 


a 33 


2£l2 





















<*12 


2e 23 










G 12 


1 









































G 23 


1 

G 13 _ 
























(37) 



Avec les conditions de symetrie 

Vi2_V2l Vl3 _^31 v 32 _ v 23 

Ej E 2 Ej E 3 E 3 E 2 

III - 2.4 Materiau elastique homogene isotrope transverse 

Un materiau homogene isotrope transverse est tel que la matrice de comportement est 
invariante par toute rotation autour d'un axe privilegie. En utilisant cette invariance, on 
montre que seuls 5 parametres independants caracterisent le comportement. Si l'axe est porte 
par la direction 3, on a alors: 





1 


-V12 


-V13 















eT 


Ei 


Ei 










Ell 




-v 2 i 


1 


-V13 

















Ei 


Ei 


Ei 












e 22 




-V31 


-V31 


1 













a 22 


e 33 




E 3 


E 3 


E 3 


1 






< 


°33 


2£i2 





















a 12 


2e 23 










G 12 


1 

G 13 






a 23 


2ei3, 


















1 

G 13_ 




°13. 






















(38) 
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III - 3 Theoreme de superposition 

Si (U, f , F) et (V, g, G) sont deux jeux de donnees engendrant respectivement des solutions 
u et v , alors ocu+|3v est solution du probleme de donnees (ocU + (3V,af + |3g,aF + |3G) . 

Ill - 4 Criteres de limite d'elasticite pour les materiaux isotropes 

Les criteres de resistance que nous allons definir representent des valeurs limites pour les 
contraintes maximales, et permettent de ce fait de garder un caractere elastique aux 
deformations. 

Ill - 4.1 Critere de Tresca 
II consiste a considerer de maniere independante les trois contraintes de cisaillement 
maximal du tricercle de Mohr. Soit en fonction des contraintes principales 

Sup|o No (39) 

III - 4.2 Critere de Von-Mises 

^((^i-^ii) 2 +(^i-^iii) 2 +(^ii-^iii) 2 ^o (40) 
ou encore 

^(( ll-°22) 2 +( a ll- 33) 2+ ( 22 - 33) 2 +6 ( 12 2 +°13 2 +O 23 2 ))^ O 



III - 5 Thermoelasticite 

Tout solide soumis a un ecart de temperature cherche a se dilater s'il le peut. S'il ne peut 
se dilater, alors il y a apparition de contraintes dites « d'origine thermiques ». Ce phenomene, 
pour des ecarts de temperature ST faibles par rapport a la « temperature de repos se traduit par 
une loi de comportement dans le cas general sous la forme: 



G(x,t)=C(x):(e(x,t)-aoT) (41) 



ou a represente le tenseur anisotrope des dilatations (en X" 1 ). 

Pour plus de rigueur, nous invitons le lecteur a se referer au polycopie du cours de 
Mecanique des Milieux Continus. 

Dans le cas ou le materiau est isotrope la loi se simplifie en : 

o= XTr(e)l + 2(^8-a(3X + 2(^)6Tl (42) 

Et la relation inverse : 

8= --Tr(o)l + ^^o + aSTl (43) 
E E 
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THEORIE DES POUTRES 

A partir de ce chapitre, on utilise les hypotheses des petites perturbations, du quasi- 
equilihre et de Velasticite lineaire isotrope. 

I - Definitions, hypotheses de Bernouilli 

I - 1 Definition d'une poutre 

On appelle poutre le solide engendre par une surface plane dont le centre de gravite decrit 
une courbe y, la surface S restant normale a cette courbe, avec: 

* La courbe y est appelee ligne moyenne ou fibre moyenne 

* La surface S est appelee section normale 

* Le rayon de courbure en tout point de y doit etre grand par rapport aux dimensions de S 

* Les dimensions de S sont negligeables devant la longueur de la courbe y 

* Les variations de forme et de dimension de S doivent etre progressives 

I - 2 Notations 

Considerons une poutre rectiligne de section droite constante So et de longueur Lo dans la 

configuration de reference. A cette configuration de reference on associe le repere orthonorme 
direct (0,e l5 e2,e 3 ) , tel que : 

* O est un point d'une section extremite de la poutre sur la fibre moyenne 

* ej est le vecteur unitaire porte par l'axe de la poutre 

* e 2 est un vecteur unitaire dans le plan des sections droites, de preference parallele a un 
axe de symetrie de S (s'il en existe un) 




On note S(Xi) la section droite d'abscisse Xj. Dans la configuration deformee, le point 
courant de la fibre moyenne deformee est x(X 1 ) , 

* le vecteur unitaire tangent a la deformee de la fibre moyenne 
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t (X j )=aF(X i )e 1 aveca= 



F(Xj)e 



* Le plan P(Xi) tangent en x(X 1 ) a x(S(X 1 )) , defini par le point x(X l ) et les vecteurs 

-l 



F(X! )e 3 etE 2 (Xj )=bF(X 1 )e 2 avecb= 



F(X 1 )e 2 



* Le vecteur unitaire E^Xj) normal a P(Xi) 

* Le vecteur unitaire E 3 (X 1 ) tel que (x(X 1 ),E 1 (X 1 ),E 2 (X 1 ),E 3 (X 1 )) soit un repere 
orthonorme direct. 




Tf 3 (Xi) 

On notera u(X)=u(X 1 ,X 2 ,X 3 ) le deplacement de la particule X, et u(Xj) celui de la 
particule situee sur la fibre moyenne. 

I - 3 Hypothese de Bernouilli 

Le caractere lineique de la geometrie des poutres fait qu'on s'attend a ce que les 
phenomenes preponderants soient essentiellement longitudinaux. On ne s'interessera done pas 
aux deformations de sections droites. On enonce alors les hypotheses de Bernouilli 

(i) Les sections droites restent planes 

(ii) Les sections droites se deforment librement dans leur plan 

(Hi) La variation des deformations de la section le long de la poutre est tres petite 

Remarques: 

* D'apres (i) on peut confondre P(X!) et x(S(Xi)) 

* Le deplacement de la section droite peut etre represente par un vecteur translation (par 
exemple u^Xj) et par un vecteur rotation (par exemple le vecteur rotation ro(Xj) du repere 
(x(X 1 ),E 1 (X 1 ),E 2 (X 1 ),E 3 (X 1 )) par rapport au repere (0,e l5 e 2 , e 3 ). 



* Le deplacement d'un point courant de la section consideree, du a la deformation de la 
section, est done de la forme 

v 2 (X)E 2 (X 1 )+v 3 (X)E 3 (X 1 ) 

les fonctions v 2 et v 3 sont nulles en Xi, et d'apres (iii) leurs derivees v 2 j et v 3 j petites 
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devant Ui , (Qj et Uij , (Oij. On peut done ecrire: 

VXe S(X X ) u(X) = u(Xj) + ©(X,) aX^X + v 2 (X) E 2 (X,) + v 3 (X) E^XJ (1) 

L'hypothese des petites perturbations fait que les composantes Ui, Vj, et C0i sont petites; ceci 
implique que E 2 (X l ) est de la forme e 2 + r\ , ainsi que E 3 (Xj ) . 

Si on explicite (1) dans la base (e 1? e 2 , e 3 ) en ne retenant que les termes d'ordre 1, on 
obtient: 





u 1 (X 1 )+0) 2 (X 1 )X 3 -co 3 (X 1 )X 2 ^ 




VXeS(X 1 )u(XH 


u 2 (X 1 )-co 1 (X 1 )X 3 +v 2 (X) 


(2) 




u 3 (X 1 )+co 1 (X 1 )X 2 +v 3 (X) j 





On remarque alors, en ne retenant que les termes d'ordre 1: 

- F(Xj ) e 2 et x(X t ) x(X 1 ,l,0) sont egaux et unitaires, done egaux a E 2 

- F(Xj ) e 3 est unitaire et orthogonal a F(Xj)e 2 
Par contre, 

f(Xj )e! )• (f(Xj )e 2 )=u 2>1 (Xj )-co 3 (X 1 ) + - 

^(Xj )e! )■ (f(X! )e 3 )=u 3>1 (Xj )+t0 2 (Xj ) + • • ■ 

Done en general la deformee d'une section droite n'est pas, au second ordre pres, 
orthogonale a la deformee de la fibre moyenne. 

II - Deplacements et forces generalises 

La geometrie des poutres fait que les sollicitations exterieures peuvent etre considerees 
comme donnees 

- soit sur une partie de la surface laterale [Xj;X} ]edS i e {0,...,J} avec X"=0 etXj=L 

- soit sur les sections extremites 

- soit sur des cercles T 1 = { X} }edS 

II - 1 Deplacement generalise 

Afin de bien distinguer le deplacement du milieu continu de celui de la fibre moyenne, on 
note 

u(X 1 ,0,0)=u f 

D'apres (2), se donner u satisfaisant a l'hypothese de Bernouilli equivaut a se donner 
u p =(u f ,0») 

On appellera u p le deplacement generalise de la poutre. 



I.S.I.T.V. 



-28- 



Theorie des poutres 



II - 2 Puissance virtuelle des efforts exterieurs 



Soit 8u une vitesse virtuelle de emplacement et § un systeme de forces exterieures, e'est a 
dire le couple ( f ,F ) densite volumique de forces et densite surfacique de forces. La puissance 
virtuelle des efforts exterieurs s'ecrit alors: 

n ext (5ui)=jn n f-5udx+n 3n F.5udx 

ou encore 



n ext (5ui)=Z y (fj s f • SudxJdXi+lY 

i=l Xj i=l x; V9S ) i=lr; 



J _ 



JF-Sudx dXi+ZjF-Sudx 



Soit Su p = (Su ,8co) une vitesse virtuelle generalisee de la poutre, et associons lui la 
vitesse virtuelle de deplacement 

5u{ + 80)2X3 - 8u) 3 X 2 
8u 2 - 8(o[X 3 



8u(X) = 80* (X,) + SoXX,) aX,X = 



8u[ 




8(o[ 








<8u 2 


> + < 


8(0 2 




x 2 


. = < 


8u f 3 




8(0* 




x 3 , 





Su* + 8(o[X 2 



et en utilisant la propriete du produit mixte ((a a b) ■ c = a ■ (b a c)), on a finalement 



j-i x i+i 



n ext (8u,^)) = 1180* 



JJ s f dx+ jFdx 

as 



dX, 



J-l X H 



+ Z j 8w 



1=1 Xj 

J 



UsXjXAf dx+ jXjXAFdx 

3S 



dX, 



(3) 



+ Z^8u t (X' 1 )- JFdx + 8uXX;)- jXjXAFdx 

i=i r r, 



II - 3 Forces generalises 



En posant 

f f (X,) = JJ s f dx+ jFdx F' = {Fdx 

as n 



c f (X 1 ) = jj s X 1 XAf dx+ JX,XAFdx C = jX^AFdx 

as n 



Soit 



n ext (8u^)=n ext (8u p ^ p )=j L 8u f (X 1 )^ f (X 1 )dX 1 +Z8u f (xj)-$ i (4) 

i=l 



avec 



y , • • • , * J Jbt^ =(f f ,c f \& =(? .C 1 ) 



On appelle (|) P force generalisee appliquee a la poutre. Cette force est constitute d'une 
densite lineique de force generalisee (]) repartie le long de la fibre moyenne et de forces 
generalisee concentrees <J>\ 
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III - Deformation et contraintes generalises 

III - 1 Deformations generalisees 

D'apres (2) en utilisant l'hypothese des petites perturbations, et en negligeant les derivees 
de v 2 et V3,on obtient 



u i,i + C0 2,iX 3 - 0) 3 jX 2 Sym 

| 3,1 +(0 1,1 X 2+^2) ^( v 2,3+ v 3,2) 



Sym 
Sym 



'3,3 



(5) 



En introduisant alors les six quantites 



ajCX^uUX!) 



X 1 (X 1 )=to u (X 1 ) 



a 2 (X 1 )=u t 2il (X 1 )-a)3(X 1 ) % 2 (X 1 )=co 2J (X 1 ) 
a 3 (X 1 Hi| >1 (X 1 ) + (D 2 (X 1 ) % 3 (X 1 )=(03,i( x i) 
on trouve de maniere simple 



(6) 



8= 



a i+%2 X 3-%3 X 2 Sym Sym 
-(a 2 -Xi X 3) v 2,2 Sym 

|(a 3 +Xi x 2 ) |( V 2,3 + V 3,2) v 



3,3 



(7) 



Les six quantites definies par (6) constituent la deformation generalisee de la poutre en la 
section S(Xi). 

Afin d'interpreter mecaniquement cette definition nous allons etudier successivement les 
cas ou une seule de ces quantites est non nulle. 



Casl: a^O, a 2 =a 3 =x 1 =x 2 =x 3 =0 

On en deduit 
u 2 f =0,u f 3 = etro = 6 
soit 

u 1 (X)=a 1 X 1 et u 3 (X)=u 2 (X)=0 

La poutre est dans un etat d'allongement 

pur 



a x L 
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Cas2: a 2 *0,a 1 =a 3 =x 1 =x 2 =x 3 =0 

On en deduit 

uf =0,113 = Oetoi> = 

u^Xj) = a 2 X 1 etu 2 (X) = a 2 X l 

La poutre est dans un etat de glissement 

dans le plan (e^ej) 

Cas3: a 3 *0, a 1 =a 2 =x 1 =x 2 =x 3 =0 

La poutre est dans un etat de glissement 
dans le plan (e l5 e 3 ) 



Cas4: y^O, a t = a 2 = a 3 = % 2 = % 3 =0 

On en deduit 

a i = Osoitu = 

et(0 2 = (0 3 = O^tu)! = x 1 X l 

La poutre est dans un etat de torsion autour 

de son axe 

Cas5: % 2 *0, a 1 =a 2 =a 3 =x 1 =x 3 =0 

uj = u' 2 = et = co 3 = 

Soit Q) 2 = x 2 X, et u f 3 (X,) = 
D'ou u 1 (X) = x 2 X 1 X 3 , 
u 2 (X) = 0,u 3 (X) = -|% 2 Xf 

La fibre moyenne se deforme selon une 
parabole dans le plan (e 1 ,e 3 ), la section 
droite tournant de % 2 ^i autour de e 2 . 

Cas6: a 1 =a 2 =a 3 =x 1 =x 2 =0 

C0j = (0 2 = et par suite u t (X) = -% 3 XjX 2 , 

u 3 (X) = 0etu 2 (X) = |x 3 Xf 

La fibre moyenne se deforme selon une 
parabole dans le plan (e 1 ,e 2 ), la section 
droite tournant de % 3 X 1 autour de e 3 
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On en deduit alors la signification des deformees generalisees: 

- a 1 est l'allongement unitaire de la fibre moyenne 

- a 2 et a 3 sont des glissements dans les plans (e 1 ,e 2 ) et (e 1 ,e 3 ) 

- Xi est Tangle de torsion autour de l'axe e 1 par unite de longueur 

- % 2 , %3 sont l es courbures de la fibre moyenne dans les plans (e 1 ,e 2 ) , (e 1 ,e 3 ) 



III - 2 Puissance Virtuelle des efforts interieurs 



En tenant compte de la cinematique particuliere des poutres, nous allons expliciter 
l'equation 

n int (o1,5)=-jnnS:o1dX (8) 

ou 8e designe une vitesse virtuelle de deformation, et a l'etat de contrainte dans la poutre. 
Nous choisirons, bien sur, 8e de la forme 



8e= 



8aj +8% 2 X 3 -8% 3 X 2 Sym Sym 

|(8a 2 -8x 1 X 3 ) 8v 22 Sym 



^{8a 3 +d% l X 2 ) |(5v 23 +5v 32 ) 8v 



3,3 



(9) 



Pour que les contraintes satisfassent l'hypothese (ii) de Bernouilli, il est necessaire que 
O n o 12 o 13 

(10) 



0= 



a 12 
°13 













En portant (9) et (10) dans (8), on trouve finalement 

-n int (s£,a) = j L [8a 1 jj s a 11 dX + 5a 2 jj s a 12 dX + 5a 3 jj s a 13 dx]dX 1 
+ !o L [SXi!! s (X 2 a 31 -X 3 o 21 )dx]dX 1 
+ ! L [S%2!lsX 3 a 11 dX-5x 3 n s X 2 a 11 dx}iX 1 

Posons alors 

T\ (Xj ) = jj s Oj jdX Mj (Xj ) = H s (X 2 o 31 - X 3 o 21 )dX 
T 2 (X 1 ) = jj s o 12 dX M 2 (X 1 ) = jj s X 3 o 11 dX 
T 3 (X 1 ) = n s o 13 dX M 3 (X 1 ) = -JJ s X 2 a 11 dX 



(11) 



(12) 



avec ces notation, on obtient plus simplement 



-n int (5E,o)=j L 



Z(5a i T i (X 1 )+6x 1 M i (X 1 )) 
i=i 



dX, 



ou encore 



T = jj s ae,dX et M = JJ S X,X a oe^X 



(13) 



I.S.I.T.V. 



-32- 



Theorie des poutres 



III - 3 Contraintes generalisees, equation d'equilibre 

III - 3.1 Contraintes generalisees 

On definit les contraintes generalisees comme etant la fonction 
(s p ) : X, e [0,L] -> (s^X,)) = (T 1 (X 1 ),T 2 (X 1 ),T 3 (X 1 ),M 1 (X 1 ),M 2 (X 1 ),M 3 (X 1 ))e M 6 (14) 

On appelle T 1 l'effort normal, T 2 et T 3 les efforts tranchants, M t le moment de 
torsion, M 2 et M 3 les moments de flexion. 

Si on considere la partie de la poutre a gauche de la section S(Xj), alors la normale unitaire 
sur S(Xj) sortante est n =6 1 ; la densite surfacique de force exercee par la partie droite de la 

poutre sur la partie gauche est ah, c'est a dire le vecteur de composante o ll ,o 12 ,o 13 . En 

consequence, 

la contrainte generalisee s p (X 1 ) est constitute des elements de 
reduction, au centre de la section S(X X ), du torseur des forces 

appliquees par la partie droite de la poutre sur la partie gauche! 

Certains etats de sollicitation elementaires sont appeles 'sollicitations simples', ils 
correspondent a des cas de chargement frequemment rencontres. Leur etude permet par le 
theoreme de superposition l'etude de cas plus complexes appeles 'sollicitations combinees'. 



Ti 


T 2ou3 




M 20 u3 


Designation 


^0 











Traction ou compression simple 





^0 








Cisaillement pur 








*0 





Torsion pure 











*0 


Flexion pure 





*0 





^0 


Flexion simple 


*0 








*0 


Flexion composee 
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III - 3.2 Equation d'equilibre 

Nous allons maintenant appliquer le principe des puissances virtuelles pour determiner les 
equations d'equilibre. Pour toute vitesse virtuelle de deplacement generalise, on doit avoir 

n ext + n int = o 

Nous avons d'apres (4) 

n ext = Z Y ( 55f ■ f f (Xj) + 8uV c f (X^JtlXj + l{8u f (Xj) ■ F 1 +8u)(X 1 1 ) • C 1 } 



i=i 



i=i 



et d'apres (13) 
nint =- Jo 



Z(5a 1 T I (X 1 )+5% i M i (X 1 )) 
i=i 



dX, 



En utilisant la definition des deformations generalisees (6), 

n int =-j L [5ui -f (X0-5(O 3 T 2 + 5(0 2 T 3 +563^ -MCXoJdX! 

Soit en remarquant que e 1 a T = T 2 e 3 - T 3 e 2 

n int =-! L [Sufi • T(Xj ) + 5©f! ■ M(Xj ) - 5co f ■ (Sj a f )JdXj 

On decompose alors l'integrale en somme d'integrales comme pour les efforts exterieurs 

nint—ZjJ^lsuf! T (XO + Srofj ■MiXJ-btf -(ej aT)^ 
i=i x i 



En effectuant une integration par partie et en notant X} la valeur de Xj pris par valeur 
superieur, 



n 



int 



= Z6 +1 [ou f • T^Xj) + 5(b f ■ M iCXO + 5(Q f ■ (ej a f )]dX 1 
i=i x i 



j-i 

- I 

i=i 



5u f ■ 



T(Xj +1 )-f(xi + )V8u) f -fM(xi +1 )-M(xj + ) 



V 



V 



Soit 



n 



mt 



= z 



J /^[su* •f 1 (X 1 ) + 6© f -M^X^ + Sa^ -(ej A^JdXi 



i=i x i 



J 

+ z 

i=l 



5u f ■ 



f(Xj )-T(Xj ) 



+ 8u) f ■ 



M(Xj )-M(xj ) 



Enfin en application du principe des puissances virtuelles 

ZjT 1 [ 5Qf ■ (f f + f j )+ 8u) f ■(c f +M 1 +e 1 A f jJdXi 



i=i x i 



+ z 

i=l 



5u f 



F 1 +T(Xj ' ) -T(Xj )] + 8a> f - C 1 +M(X}' )-M(X} ) 



(15) 



= 



Com me (15) est verifie pour toute vitesse virtuelle, on obtient pour l'equilibre 



Pourchaqueint ervalle X \ , X J +1 
PourXji = l,...,J 



f f +f j 







(a) 

c 1 +Mj +ej aT =0 (b) 
F i +f(X 1 1 + )-f(X 1 1 ") = (c) 
C i +M(X 1 1 + )-M(Xi") = (d) 



(16) 
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Remarque importante: 

Sur un troncon de poutre JxJ , X J" 

f j =6etM 1 +e 1 aT = 6 
soit en particulier 

dM2 = T 3 et^=-T, 



i+i 



dans le cas ou il n'est pas charge, on a 



dX 



l 



dX 



IV - Loi de Comportement elastique lineaire 



D'apres (10), a cause de l'hypothese (ii) de Bernouilli a= 



Ojj Sym Sym 
o 12 Sym 
G 13 

Or, pour un materiau elastique la loi de comportement en fonction du module d'Young E et 

du coefficient de poisson est 
= v = = l+v = 



Soit dans notre cas 
1 



'11 - ~°11' £ 22 ~ £ 33 
E 



-a,,,8 



U' c 12 



1+V _ i+v 

~~ ~°12' E 13 - ~ ^°13' e 23 _ (J 



D'ou 



T 1 =n s a 11 dX = n s E8 11 dX = n s E(a 1 +x 2 X 3 -X3X 2 )dX 
T 2 =JJ s a 12 dX = n s ^8 12 dX = n s ^-|-^(a 2 -%!X 3 )dX 

T 3 = jj s a 1 3dX = n s ^-E 13 dX = jj s ^|^(a 3 +%1 X 2 >iX 

Mj =jj s (X 2 o 31 -X 3 o 21 )dX = jj s — ^-(-a 2 X 3 +a 3 X 2 + %1 (X 2 ! +X 2 3 ))iX 

M 2 = n s X 3 o n dX = jj s E( ai X 3 + %2 X 3 2 - %3 X 2 X 3 )dX 

M 3 =-jj s X 2 o 11 dX = -jj s E(a 1 X 2 - %3 x2 +%2 X 2 X 3 )dX 

Pour simplifier les ecritures, nous definirons la fibre moyenne telle que 
JJ S E(X)OMdX=0 

c'est a dire 

JJ s EX 2 dX = JJ s EX 3 dX=0 

de plus nous choisirons e 2 de telle maniere que 
j]" s EX 2 X 3 dX = 

Dans ces conditions, nous obtenons 
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T 2 =a 2 jJs— ^— dX-xJJs^^X.dX 



2(1 + v) 



■2(l + v) 



2(1 + v) 



Js 2(l + v)^ 2 



Ml=j!s 2(i+^) %l(X ' +X ' )dX 

M 2 = %2 {j s EX 3 2 dX 
M 3 =%3!!sEX^dX 

Si de plus les caracteristiques ne dependent pas de l'espace, c'est a dire que Ton peut sortir 
E et v d es integrates, on obtient: 





Tl =a lE S,T 2 = ^ ES ,T 3 = a3ES 
2(1 + v) 3 2(1 + v) 


(17) 


et 


M 1 =-^-jj s (x2 + x2)dX,M 2 =E %2 jj s X 3 2 dX,M 3 =E %3 n s x2dX 




ou encore 






M 1 =^^y,M 2 =E %2 I 2 ,M 3 =E %3 I 3 


(18) 



ou Ij est le moment d'inertie autour de l'axe (O,^) . 
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ETUDE DE SOLLICITATIONS SIMPLES 



I - Traction ou compression 



I - 1 Definition 

On dit qu'une poutre est dans un etat de traction (ou compression) quant le torseur des 
actions exterieures est de la forme: 

T^O, T 2 = 0, T 3 = 0, M 1= 0, M 2 =0, M 3 =0 

f f =f f e l5 c f =6,^ = F i e 1 ,C i =0 
Attention : lorsque la longueur est superieure a environs 8 fois la plus grande dimension 
transversale, une poutre sollicitee en compression est calculee au "flambement". 

I - 2 Deformations et contraintes 



De (11-17) et (11-18) on deduit: 

T^ajES, a 2 = a 3 =0et %i = Xi = %3 = 
Le tenseur des contraintes qui satisfait l'equilibre est de la forme: 



a n 



(1) 



(2) 






ou G n est une valeur constante dans toute la section. On obtient alors pour le tenseur des 
deformations: 



8= 



—2- 
E 

-v^ 
E 















-v^- 
E 



(3) 



Relation de la contrainte avec 1' effort normal 

On sait que 
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Ti 

cr n=— 
S 



(4) 



Allongement de la poutre 

Soit AL l'allongement subi par une poutre de longueur L. Par definition 
AL=f Y d Ul 



Dans le cas present, le deplacement d'une section droite est une translation d'axe e 1 . 




Si la section de la poutre est constante 



AL= T > L 




(7) 


SE 



II - Torsion 

Les sections droites de contour quelconque, lorsqu'elles sont sollicitees en torsion, se 
gauchissent. Ce phenomene remet en cause l'hypothese de Bernouilli. Les poutres droites de 
section circulaire ne subissent pas de gauchissement. 



II - 1 Definition 



Une poutre est sollicitee en torsion pure si: 

T 1 = 0, T 2 = 0, T 3 = 0, M^O, M 2 =0, M 3 =0 (8) 

f f =6, c f =f f e„ F=6, C=C% 



II - 2 Deplacement, contraintes, deformations 

De (8), (11-17) et (11-18) on deduit: 

aj = 0,a 2 = 0,a 3 = Q, Xl ± 0,% 2 = 0,% 3 = (9) 

D'apres (9), on a u = , la section tourne done uniformement autour de son axe. %i est 

l'angle unitaire de torsion. 

Comme 8 n = aj +% 2 X 3 -% 3 X 2 = alors o n =0 , soit a= 



o 12 o 13 
o 12 
o,. 
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De la loi de comportement, £ y - = ^r^y _ g'^kk^ij > on ti re 




l+v 



13 










X1X3 

2 

XiX 2 



X1X3 XiX 2 



2 2 




Dans le repere(G,e 1 ,e r ,e e ) 
V M u(M) = -0) 1 X 3 e 2 + 0) iX 2 e 3 
Or e r = CosG e 2 + Sin0 e 3 
et e e = - SinG e 2 + CosG e 3 

Soit u(M)=0) 1 re e 



(10) 



Pour determiner les deformations, nous utilisons la definition des deformations en 
coordonnees cylindriques: 



£(G,ei,e r ,e e )- 



3uj 

Sym 
Sym 



1 (d"i I 

2 V 3r ^ 3X, 

Sym 



) 1 (l d"i 1 3u e ) 

/ 2 Vr 36 3X, / 



+ 



9ug 
3r 






Sym Sym 



3u P 



1 

2 vax 



Soit, en utilisant (II-6) 











2 










rXi 

2 


£(G,e 1; e r ,ee) _ 
























Sym 










Sym 









et 



a(G,ei,e r ,e„) = 



X 



TOO 



avec T = 2Ge Xi e et G = ^^ 



G est appele le module de glissement ou module de Coulomb. En definitive, on peut ecrire: 
1 = 0% (11) 



Relation entre la contrainte et les efforts generalises 

En utilisant (11-18) et (11) M { 



Soit 




2(1 + v) r 



Mi 
x = — -r 
h 



(11) 



La contrainte de cisaillement est maximale a la peripheric. 
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Determination de Tangle de torsion 

Xi est Tangle unitaire de torsion. Done, si on note AB Tangle de torsion entre deux sections 
A et B, on obtient: 

e AB =g A B Xidx 1= gB JL^ 

soit 



0AB = f x B ML dXi 

Xa GI i (12) 



II - 3 Exemple 

Un arbre en acier de longueur L=lm est sollicite en torsion par un couple M=1500 mN. 
Sous V action de ce couple, on desire que V angle unitaire de torsion reste inferieur a une 
valeur limite OCi=0.25 °/m et que la contrainte de cisaillement soit inferieure a 120 N/mm 2 . 
On prendra G=8.104 N/mm 2 et p=7800 kg/ml 

a) Calculer le diametre admissible Dj de I 'arbre. 

b) On suppose que V arbre est un tube de diametre exterieur D e =90 mm. Quel doit-etre le 
diametre interieur ? 

c) Quelle economie de masse a-t-on realisee ? 
Reponses: 

a) Pour que la geometrie de Tarbre soit admissible, il est necessaire de satisfaire deux 
cri teres. 

* Critere de deformation 

Onveut %i < ocl soit %, < a L ou encore L> =43.10 mm 

GIj Ga L 

Or I 1= ff (x 2 2 + X 3 2 )dX=nr 2 rdrd0=2;r^=7C^- 



Soit Dj> 



32M 



7cGa L; 



4 

=82mm 



* Critere de resistance 



On veut T=^^-<120N / mm 2 soitD 1 > 



2Ii 



16M 
12071 



3 



~40mm 



Pour satisfaire au cahier des charges, Tarbre doit au moins avoir 82 mm de diametre. 

b) Pour un tube, le moment quadratique est l l = 7C — . 

Avec la condition I 1 <43. 10 5 mm 4 , on trouve Dp68.3 mm 

c) II est immediat de constater que le gain de poids est de environs 20 Kg. 
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III - Flexion 

III - 1 Flexion pure 

Une poutre est sollicitee en flexion pure lorsque: 
Ti=T 2 =T 3 =M 1 =0 et M 2 *0 , M 3 *0 

De (11-17) et (11-18) on deduit: 

a 1= a 2 = a 3 = et %, =0, % 2 *0, % 3 * 0, 
SoitE n = % 2 X 3 -% 3 X 2 ete 12 = £ 13 = 

a n o 12 a 13 



Comme = -fa^eta = 



a 12 
°13 








On obtient: 

a 12 = a 13 = Oeta u = E% 2 X 3 - E% 3 X 2 

Soit, d'apres (11-18): 



(13) 



n M 2 






<^ll - "J - X 3 


I 

A 3 


(14) 



III - 2 Flexion pure plane 

Une poutre est sollicitee en flexion pure plane lorsque: 
T 1 =T 2 =T 3 =M 1 =M 2 =0 et M 3 *0 

D'apres les developpements du paragraphe precedent: 



o n =- 



I, 



"X 2 - -E% 2 X 3 



(15) 



La contrainte est maximale a la peripheric de la poutre. 



Zone comprimee a<0 



Zone tendue a>0 
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Deformee de la poutre 

D'apres (II-6), on peut ecrire: 
a 1= -> ufi=0 
% 1= » u =0 

% 2 =0 »2,1 =0 
a 3 =0 u^=0 
a 2 =0 u 2a =(0 3 

Soit, en utilisant (11-18) 



uf =Cste = 
— > C0j =Cste = 
-> (0 2 = Cste = 
u 3 = Cste = 



"2,11 = (0 3,1 =%3 



2..f 



d z u 



dX; EI 3 



(16) 



III - 3 Flexion plane simple 

II s'agit du cas particulier oil: 

T 1 =T 3 =M 1 =M 2 =0 et M 3 *0 , T 2 *0 

(ou Ti=T 2 =Mi=M 3 =0 et M 2 *0 , T 3 *0) 
La poutre n'est soumise qu'a des efforts tranchants et des moments flechissants. 



(17) 



* Contrainte normale due au moment flechissant 

°n = E( ai + % 2 X 3 - E% 3 X 2 ) = -E% 3 X 2 , soit: 

M 3 



I- 



X. 



(18) 



* Contrainte tangentielle 

A 




Xi Xi+dXi 

Considerons un troncon de poutre, non-charge, compris entre Xi et Xi+dXi. 



Le probleme est plan de telle sorte que a = 



a, 1 o 12 



11 

°12 










Les equation d'equilibre nous donnent a lu + o l2 2 = , o 12 x = et T j = 

T u = 0h-»JJo 11>1 dX 2 dX 3 =0 
s 

JJo lu dX 2 dX3 + JJo lu dX 2 dX 3 =0 

SI S2 
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Jf-o 12i2 dX 2 dX 3 + JJ 

SI S2 



■M, 



V i 3 



X 



dX 2 dX 3 =0 



-bo 12 + JJ^-X 2 dX 2 dX 3 =0 
Soit en definitive, 

G 12 =^tfX 2 dX 2 dX 3 (19) 

bl 3 S2 

* Deformee de la fibre moyenne 

1 y 

Rappel: Pour une courbe d'equation y=f(x), la courbure est definie par — = 



R (l + y') 3/2 

Done pour des deformees "petites", e'est-a-dire y'«l, on obtient 1/R=y". 
Soit dans notre cas: 
d 2 u 2 _ M 3 



dX 2 EI 3 



(20) 



III - 4 Exemples 



III - 4.1 Exemple de flexion pure plane 
F 1 F 



A 

A" 



D 



B 



Une poutre droite rectiligne de section constante repose sans 
frottement sur 2 appuis simples en A et B et supporte une 
charge concentree F=1500 N en C et D. 



AC=CD=DB=a=0,5 m. 
Ra If If Rb| 



D 



B 



a) On cherche a determiner les reactions aux points ou sont 
imposees des conditions cinematiques (en A et B). Comme il 
s'agit d'appuis simples, il ne peut y avoir que des reactions 
(pas de couples). 

Pour determiner les reactions, on applique le principe fondamental de la statique. 
Soit Rax, Ray, Rbx, Rby les composantes des reactions en A et B. 
E Forces = OsoitR A + 2F + R B =0 



R AX 




' 0" 




(T 




R BX 







' Ray 


> + < 


-F 


> + < 


-F 


> + < 


Rby 


> = < 


0> 












j 













E Momentsen A = Osoit AC aF + ADaF + ABaR b =0 



a^ 




0^ 




2a^ 




' 




V 




R BX 







■ 




-F 


> + < 





■ A< 


-F 


> + < 





■ A< 


Rby 


> = < 


0> 












0^ 




j 




0^ 













Les equations d'equilibre et la symetrie du probleme impliquent: R A = R B =-F 
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b) On cherche les contraintes generalisees dans la poutre. 

On sait que les contraintes generalisees dans une section donnee sont egales au torseur des 
actions exterieures a droite de la section, ou encore, a l'inverse du torseur des actions 
exterieures a gauche. 

* soit une section G dans le troncon de poutre AC. 

-f = R A soit T, = T 3 = et T 2 = -F 

-M = GAaR a soitM 1 = M 2 =OetM 3 = xF 

Le tron£on AC est dans un etat de flexion plane simple. 

* soit une section G dans le troncon de poutre CD. 

— T = R A +FsoitT\ = T 2 = T 3 = 

-M = GAaR a +GC AFsoitMj =M 2 =0etM 3 = aF 
Le tron£on CD est dans un etat de flexion pure. 



Ill - 4.2 Exemple de flexion plane simple 
P=1500N P=1500N 



C 
lm 



A 

7T 



B 



D 



2m 
Ra P 



4 



'I 



2m 
Rb 



C A 



i 



D 



B 



Une poutre droite rectiligne de section circulaire 
constante repose sans frottement sur 2 appuis simples en 
A et B. La poutre est constitute d'un materiau de limite 
elastique de 1600 bars. Quel doit etre le rayon de cette 
poutre ? 

a) On cherche a determiner les reactions aux points ou 
sont imposees des conditions cinematiques (en A et B). 
Comme il s'agit d'appuis simples, il ne peut y avoir que des 



reactions (pas de couples). 

I Forces = OsoitR A + 2P + R B =6 
R A +R B =2P 

X Momentsen A = Osoit AC a P + AD a P + AB a R B =0 
P-2P + 4R B =0 
Soit R B =P/4 et R A =7P/4 

b) On cherche les contraintes generalisees dans la poutre. 

On sait que les contraintes generalisees dans une section donnee sont egales au torseur des 
actions exterieures a droite de la section , ou encore, a l'inverse du torseur des actions 
exterieures a gauche. 

* soit une section G dans le tron§on de poutre CA. 

- M = GC a PsoitMj = M 2 = 0etM 3 = -xP 

* soit une section G dans le tron§on de poutre AD. 
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— M = GC a P + GA a R A soitMj =M 2 =0etM 3 =-Px + 7P(x-l)/4 

* soit une section G dans le troncon de poutre DB. 

M = GB a R B soitM 1 = M 2 = OetM 3 = P(5- x)/4 

Le moment maximal est atteint pour x=l, soit M 3max =-l*P. 

4M, 





r M 3 ^ 




f 4M 3 N 




Max 




= Max 


T X 2 


d'ou R = 3 - 








VTlR 4 J 


\ 



no. 



III - 5 Etude de la deformation des poutres en flexion 



III - 5.1 Methode de la double integration 



On a vu que 



U f 2=! 



d\£ 

dXj 2 



M : 
EL 



done par integrations successives on obtient: 



M 5 1 
J-^dXi dX 1 +C 1 X 1 +C 2 

. EI 3 J 



(21) 



ou Ci et C 2 sont des constantes qui seront determinees par les conditions aux limites. 
Exemple: 



A |b 

A 



F=16kN 

C 



D 



Ra 



F 



Rd 



D 



lm lm 2m ^ 

On cherche a determiner la deformee de la poutre. 

a) Premierement on cherche les reactions aux points d'appuis. 

Jr a +R d +F = 6 |R a +R d =F|R a =3F/4 
|aBaF+ADaR d =0[-F+4R d =0[R d =F/4 

b) On determine maintenant le moment de flexion dans chacun des tron§ons de la poutre. 
* Dans le tron§on AB: 0<x<l 

M = — GA a R A soitM 3 = R A X = 3FX/4 = 12X 
Soit ElV = 2X 3 +C,X + C, 



* Dans le tron5on BC: l<x<2 
A = -GA a R A - GB a FsoitI 

Soit EI 3 u 2 =-2(X-4) 3 /3 + C 3 X + C 4 



M = — GA a R A — GB a FsoitM 3 = R A X + F(1-X) = 16-4X 



c) Pour determiner les constantes on utilise les conditions aux limites et les conditions de 
continuite. 
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u' 2 



AB 



AB 



,(0) = 
(l) = U f 2 



BD 



(1) 



l 2,l 



AB 



(D=4i 



BD 



(1) 



cequidonne< 



4 



BD 



(4) = 



Cj =-14 
C 2 =0 
C 3 =10 
C 4 =-40 



En definitive u f 2 (X) = 2X 3 -14Xetu^ 

/ ab L 



BD 



(X) = -2(X-4r +10X-40 



III - 5.2 Fonctions de singularity 

Les fonctions de singularity permettent d'exprimer analytiquement une discontinuity. 
On definit la fonction de singularity d'ordre n: 

Sin < f n (X) = oolorsqueX = a 
f n (X) = OlorsqueX * a 
Sin > f n (X) = (X - a) n lorsqueX > a 
f n (X) = OlorsqueX < a 



f n (X) = (X-a) telleque 



De la meme maniere on definit les regies d'integration suivantes: 

X ' .n-l 



lorsque n<0 |(x-a) dx = (X-a 



^ \n ( X ~ a ) 

lorsque n>0 j (x - a) dx = - 

n + 1 



n+l 



Utilisation pour le calcul des fleches 

Soit q(X) un chargement vertical agissant sur une poutre. En accord avec les notations 

utilisees dans ce cours (11-16), on a: 

T 2 =-|qdX, et d'apres (II-16b) M 3 =-|T 2 dX puis en utilisant (16), on determine la 

fleche. Done, par quatre integrations successives on determine la deformee de la poutre. 
Principales fonctions de singularites et leur utilisation 

\ a 



V q = -M (X-a 



Wo 



q = W (X-a 
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Utilisation dans l'exemple precedent: 

On sait que R A =12kN et R B =4kN. 

D'ou q(X) = 12(X) _1 - 16(X - 1)" 1 + 4(X - 4)" 1 

Par integration, T 2 = -12(x)° + 16(X - 1)° - 4(X - 4)° + cste 

La constante est nulle car pour X<0 et X>4 il n'y a aucun effort tranchant. 
Par une nouvelle integration M 3 = 12(x) ! - 16(X - 1) 1 + 4(X - 4) 1 + cste 

La constante est nulle car pour X<0 et X>4 il n'y a aucun effort moment flechissant. 
Et par la suite: 

EI 3 U2,ii =M 3 =12(X> 1 -16(X-l> 1 +4(X-4) 1 

EI 3 u 2 ,1 = 6(X> 2 - 8(X - 1> 2 + 2(X - 4) 2 + Cj 

EI 3 u 2 =2(X) 3 -8(X-l) 3 /3 + 2(X-4) 3 /3 + C 1 X + C 2 

Pour determiner les deux constantes, on utilise les conditions aux limites; le deplacement 
vertical est nul pour X=0 et X=4. Soit Ci=-14kN.m 2 et C 2 =0. 



Ill - 5.3 Poutres constituant un systeme hyperstatique 

Jusqu'alors, nous n'avons etudie que des poutres formant des systemes isostatiques, ou 
statiquement determines; c'est-a-dire que nous pouvions determiner les reactions a l'aide des 
seules equations du principe fondamental de la statique. Les reactions qui ne peuvent etre 
calculees par les seules equations d'equilibre determine le degre d'hyperstaticite d'un systeme. 

I P Nous avons dans ce cas, 2 equations d'equilibre et 2 
W inconnues Ra et Ma: le systeme est isostatique. 

B 
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Ma 



Ma 




Nous avons dans ce cas, 2 equations d'equilibre et 3 
inconnues R A , R B et M A : le systeme est hyperstatique 

de degre 1 . 

Nous avons dans ce cas, 2 equations 
d'equilibre et 5 inconnues R A , R B , Rc 
Mc et M A , M c : le systeme est hyper- 
statique de degre 3. 



Exemple: Resolution en utilisant les conditions geometriques 



Rc 




\ A B C 



= Ma 



* Conditions d'equilibre statique: 

R A +R c -P=0 

M A +(a + b)R c -aP = 

* Equation de la deformee 

q(X) = -M A (X)" 2 + R A (X)" 1 - P(X - a)" 1 + R c (X - (a + b))" 1 

T 2 (X) = M A (X)" 1 - R A (X)° + P(X - a)° - R c (X - (a + b)>° 

M 3 (X) = -M A (X>° +R A (X) 1 -P(X-a) 1 +R c (X-(a + b)> 1 

EI 3 U^ 1 (X) = -M A (X) 1 +R A (x) 2 /2-P(X-a) 2 /2 + R c (X-(a + b)) 2 /2 + C 1 

EI 3 u^(X) = -M A (x) 2 /2 + R A (x) 3 /6-P(X-a) 3 /6 + R c (X-(a + b)) 3 /6 + C 1 X + C 2 



* Utilisation des conditions aux limites. 
Onsaitqueu 2 (0) = u 2 (a + b) = et u 21 (0) = 



Done Ci=C 2 =0 
L 2 - 

2L 2 2L 
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Exemple: Resolution en utilisant la methode de superposition 

IP JP 
a W b 



B C 




ABC 

Probleme I 



Probleme I: 

L'equilibre nous donne Rai=P et MAi=Pa 

q I (X) = -M AI (x)- 2 +R AI (x)- 1 -P(X-a)- 1 

qi (X) = p(- a(x)" 2 + (X)" 1 - (X - a v " 1 



A 

Probleme II 



Rc 



T 2I (X) = PaX 



X) + (X-a^° 



M 3I (X) = p(- a(x)° + (X) 1 - (X - a) 1 ) 



EI 3 u 21 1 (X) = Pi- a(x) 1 + (X) 2 / 2 - (X - a) 2 / 2 + C 



EI 3 u 2I (X) = p(- a(X> 2 / 2 + (X) 3 / 6 - (X - a} 3 /6 + CjX + C 2 ) 

f f 

Et en utilisant les conditions aux limites : u 2I (0) = 0etu 2I j(0) = 

EI 3 u 2I (X) = p(- a(x) 2 / 2 + (X> 3 / 6 - (X - a) 3 / 6) 

Fleche en L=a+b 

EI 3 u 2I (L) = p(-aL 2 /2 + L 3 /6-b 3 /6)=Pa 2 (-3L + a)/6 



Probleme II: 

II suffit de remplacer P par -R C i et a par L. 

EI 3 u f 2II (X) = -R c (- L(X) 2 / 2 + (X) 3 / 6) 

Fleche en L=a+b 

EI 3 u^ n (L) = R c L 3 /3 

Superposition des problemes : 
Comme il y a un appui simple en C, on doit ecrire, 
u f 2I (L) + u f 2II (L) = 

Soit 

2R C L 3 +Pa 2 (-3L + a) = et R c =^-(3L-a) 
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METHODES ENERGETIQUES 



I - Theoremes de 1' energie en elasticite lineaire 

I - 1 Notations et definitions 

On note u, £ et G les solutions du probleme general d'elasticite lineaire. Les autres champs 
de deplacement seront notes v , ceux de deformations e et ceux de contraintes s . 

a) Champ de emplacements cinematiquement admissible (C.A.): 

Un champ de deplacements v est dit cinematiquement admissible, si il satisfait: 

- les conditions de regularite (continuite et differentiabilite) 

- les conditions aux bords v(X) = U(X)VXg3Q u 

b) Champ de contraintes statiquement admissible (S.A.): 

Un champ de contraintes s est dit statiquement admissible, si il satisfait aux equations 



d'equilibre: 




div(s)+f(X)=OVXeQ 
s-n=F(X)VXe3Q F 



(1) 



c) Energie de deformation elastique: 

On appelle energie de deformation elastique d'un champ de deformations e : 
W(e)=JU n CQ(e)dXouCQ(e) = -Ce : e 



(2) 



d) Energie complementaire elastique: 

On appelle energie complementaire elastique d'un champ de contraintes s : 



W*(s)={{{ n w*(s)dXou(Q*(s) = -C s:s 



(3) 



e) Energies potentielles: 

On appelle energie potentielle elastique d'un champ de deplacements v C.A.: 



^(v)= W(Bv) -Hl n f (X) ■ v(X)dX - Jjg np F(X) ■ v(X)dX 



(4) 



ou Bv est le champ de deformations du au champ de deplacement v , e'est-a-dire: 




On appelle energie potentielle elastique d'un champ de contraintes s S.A.: 




(5) 
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I - 2 Theoreme fondamental 



Proposition: 



Pour tout champ de deformations e et tout champ de contraintes s , on a: 
W(e) + W*(s)-JU n e(X):s(X)dX>0 (6) 

et l'egalite n'a lieu que si et seulement si e et s satisfont la loi de comportement 

s = C:e. 
Demonstration: 

C est un operateur defini positif. Done, 
VX e W, y~s(X) e R 3 ® s R 3 = Max \?(X) : ~s(X) - -?(X) : C(X)7'(X)\ = -s(X) : C (X)s(X) 

e'(X)ER 3 ® s R 3 [ 2 J 2 

Comme le maximum est atteint pour e(X) = C (X)s(X), on peut ecrire 

\fe{X) e R 3 ® s R 3 , y~s(X) e R 3 ® s R 3 \~e(X) : C(X)e(X) + fs(X) : C (X)~s(X) - e(X) : ~s(X) > 
Par integration, on retrouve done bien (6). 

II est evident que si la loi de comportement est satisfaite pour e(X) et s(X), alors nous avons 
egalite. 

Theoreme fondamental: 



Le triplet (u,£,a) est solution du probleme d'elasticite lineaire 
uestC.A. 

VvCA. 



aestS.A. 

wv rC \( a ~~ i ,T-)sietseulementsi 
VX g LIe = -jr ^gradu + grad u f 



o = Ce 

Demonstration: 



VsSA. (7) 
^(v)>^(u) = ^(5)>^(s) 



Soient v C.A. et s S.A.; d'apres (6) on a 
W(Bv) + W*(~s )-j\j a Bv :~s dX >0 

Par application du principe des puissances virtuelles et des equations d'equilibre (1-18), on a 
ti\ n Bv:ldX =HlJ.vdX + g aF F.vdX + g au = sn-vdX 

= \\\J-vdX + \l aF F-vdX+\\ dav = sn-UdX 

En reportant dans (6) 

W(Bv)-\j\J-vdX-jj snF F-vdX > -W*( = s) + \\ dclu = sn-UdX 

Soit, en utilisant (4) et (5) 
Z,(v)>Z,*( = s) 
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ll - Energie de deformation en RDM 

II - 1 Cas general 

Nous sommes dans le cas ou £(u) = £*(o) ,donc : 
W = i{{{ a o : EdX = Ifff f (X) ■ u(X) dX + I JJ^ F(X) ■ u(X) dX + 1 JJ^ (o(X) ■ n(X)) ■ U(X) dX 

Soit le cas de l'elasticite lineaire isotrope 

w = ittL(°xx£xx +a yy £ yy +o zz £ zz +2o xy e xy +2o xz e xz +2o yz e yz )dX 

W = JJil [2E (°xx + °yy + °zz )- t (oxyOxy + °xz°xz + °yz°yz )+ 2G (°xy + ^xz + °yz )Jd X 



Ev 



(e 



2(i+v)(i-2v) V xx yy zz 



) 2 + g(, 



p 2 +p 2 + e 2 

°xx yy zz 



)+ 2g(i 



e 2 , + e 2 +e 2 



xy ^xz 



yz 



dX 



et pour les poutres 



ZlaiTi+XiMj) 

i=l 



dX, 



(8) 



Mais comme d'apres (11-17) et )II-18) 

l\=a 1 ES,l 2 =^ 2 ^-,l^=^^_,M 1 =^ I L,M ,=Ex,I„M,=Ex,I, 
' 2(1 + v) 2(l + v) 1 2(1 + v) - 

on obtient 



ES 



a 2 

2 "2 

af + - 



■ + - 



a 2 ^ 
a 3 



V 



2(1 + v) 2(1 + v) 



+ E 



J12CL + I % 2 +I 2 

2(1 + v) 2 2 3 3 



dX, 



(9) 



ou encore 




- + ^ + - 



- + - 



Mf M: M 



■ + ■ 



■ + ■ 



ES GS GS GI, EL EL 



dX 1 = Z(F 1 -u 1 + C 1 -9 1 ) 



(9) 



Dans le second membre de (9), pour l'exemple et par soucis de simplification, nous ne 
presentons que des chargements ponctuels, mais la generalisation ne pose aucun probleme. 

II - 2 Cas particulier de la Traction/Compression 



Soit un barreau de longueur L, de section S, soumis a une force F a chaque extremite. 
On a done, T,=F et T2=T 3 =M 1 =M2=M 3 =0 et 

W=-J L ^dX! 
2 ES 

Soit dans notre cas 



W= 



F 2 L 
2ES 



(10) 
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II - 3 Cas particulier de la flexion plane simple 

Soit une poutre droite soumise a des charges transversales, c'est-a-dire T 1 =T 3 =Mi=M2=0 
et T 2 *0 , M 3 *0 
Dans ce cas 



w=k L 



T 2 + M 2 



GS EL 



dXj 



(11) 



On peut verifier que, dans la plupart des cas, l'energie associee aux contraintes de 
cisaillement est negligeable comparativement a l'energie associee aux contraintes normales. 
Soit, 



J0 2EI, 



Exemple: 



(12) 



A 



A 



^-f( X )%p(x-i; 



M 3 = i (x) 1 -Px- 2 



jfTldXi = Io /2 ^T dx i + fc2 Z r dx i 



P 2 L 



jjM 2 dx 1= r p:x2 '- ;L p 



P 2 L 3 



dX+f, /9 — (X-L) z dX = 
4 JL/2 4 48 

Soit une section circulaire de rayon R: S = 7TR 2 I 3 =tiR 4 /4 



f^dX 



o GS U ^! 2(1+ v) 1 P 2 L EtiR 4 48 N 

x — -x x — — x - - =6(1 + v) 



10 EU 



P 2 L 3 



f-T 



«1 



Car R«L. Done on peut legitimement neglige l'influence du cisaillement. 



II - 4 Cas particulier de la torsion 



Soit une poutre droite soumise a un moment de torsion, c'est-a-dire T2=T 3 =Ti=M2=M 3 =0 
et Mi#) 
Dans ce cas 



W = f 

Jt 



l M: 



2GI, 



(13) 
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III - Theoreme de reciprocite de Maxwell-Betti 

On considere deux etats d'equilibre d'un meme systeme. 

- L'etat d'equilibre I, defini par les champs o , £ ,u T , correspond a des forces volumiques 
f 1 et a des efforts surfaciques F 1 sur dQ.. 

divo' + !'=() dans a 

o I -n=F I sur dQ. 

=n =ii _ n 

- L'etat d'equilibre II, defini par les champs o , 8 ,u , correspond a des forces 
volumiques f 11 et a des efforts surfaciques F 11 sur dCi. 

div^" + f 11 = dans Q. 

i n .n = F u sur dQ. 

Par application du principe des puissances virtuelles, au champ de contraintes o en 

prenant pour champ de deplacement virtuel u , on obtient: 

f f f a ! :l n da = f f f f 1 ■ u 11 da + ff F 1 ■ u 11 ddCl 
jjjci JjJn JJdn 

= 1 Zll =IZZ=II =11 =1 

Or o :e = e Ce =o:e,soit 

tf\J 11 ■ Q 1 da + Jj^F 11 ■ Q 1 daa = JJ^f 1 ■ a 11 da + JJ^f 1 ■ u 11 daa (14) 

Le travail d'un systeme de forces (f 1 ,^ 1 ^ dans le deplacement produit par le systeme 
de forces (f n ,F n ) est egal au travail du systeme de forces (f n ,F n ) dans le deplacement 
produit par le systeme de forces (f 1 , F 1 J. 



Application; 



Dans le probleme suivant, la fleche est donnee par 

y = — ^-(2L 3 -3L 2 X+X 3 ) 
3 6EI V } 



Quelle est la fleche, a l'extremite de la poutre si on applique une 
charge Q a une distance 'a' de cette extremite? 



Soit (yp)Q la fleche au point d'application de la charge P due a la charge Q. 
Soit (yQ)p la fleche au point d'application de la charge Q due a la charge P. 
Par application du theoreme de reciprocite: 

(yp)q p = (yq)p Q 

Soit 

(y P )o = -— ( 2l3 - 3L 2 a + a 3 ) 
Jp Q 6EI V ' 
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IV - Theoreme de Castigliano et applications 

IV - 1 Theoreme de Castigliano 

Considerons un systeme de forces P ; (i=l...n) applique a une structure. Ces forces et les 
reactions constituent le systeme I. A chaque point d'application de ces forces, on a un 
deplacement (UjX- D'apres (9), l'energie de deformation est alors 

On cherche la variation de 1 'energie par rapport a une force donnee ou encore la derivee 
de l'energie par rapport a une force par exemple P r . Par definition de la derivee, on a pour 
tout accroissement 8P r : 



5P r = lim-(W(P r + 08P r ) - W(P r )) 



3P r r ^oe' 

Augmentons done la valeur de P r , d'une quantite 08P r . On constitue ainsi un systeme de 
forces II. En chaque point d'application des forces on aura alors un deplacement u ; + 08u i , et 

w n = ^-i(P i + e5P 1 )-(u 1 + e8u i ) 

L i 

Par application du theoreme de reciprocite, on peut alors ecrire 

l^+esp^o^Lp^+eSu,) 

i i 

Soit 

6P r u r = IP 1 -8u i 

i 

Puis 

Wn-w, =^-Z(P i +e5P i )-(u i +e5u i )-izP i -fl i 

Z i Z i 

= ^ZP i -05u i + ^Z05P i -(u i + e5u i ) 

2 i Z i 

= -08P ■ u r + -05P ■ (u r + 08u r ) 

2 2 

= -0 2 8P -8U+08P -u r 
2 

Ou encore 



dP e^o 



-0 2 8P -8u +08P -fl =8P Q 



Done avec un abus de notation, on ecrira classiquement: 

^r=u r (is) 

ou u r est la valeur du deplacement compte positivement dans le sens d'application de 
la force 
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En generalisant le developpement precedent, on peut enoncer le theoreme de Castigliano: 
La projection du emplacement du point d' application d'une force sur la direction 
de cette force est egale a la derivee partielle de l'energie de deformation par 
rapport a cette force. 

Le vecteur rotation du point d'application d'un couple quelconque, projete sur 
l'axe de ce couple, est egal a la derivee partielle, par rapport au moment de ce 
couple, de l'energie de deformation. 



IV - 2 Consequence: Principe du travail minimum ou theoreme de 
Menabrea 



Considerons une poutre hyperstatique reposant sur des appuis invariables. 

R'2i 

N Mb 




Ma 



Les appuis introduisent 6 inconnues R l5 R 2 , M A , R\, R' 2 , M B . Or, il n'y a que 3 equations 
d'equilibre; le systeme est done 3 fois hyperstatique. 

Rendons la poutre isostatique en supprimant les liaisons surabondantes (par exemple en 
A). 



Mb 




Ma 

Ce systeme a la meme energie de deformation que precedemment. Par application du 
theoreme de Castigliano, comme la section A est encastree, on obtient: 

~ = , — — =o , =0 

dR, dR 2 3M A 

ce qui donne 3 equations lineaires en Ri, R 2 , Ma- 



Les valeurs que prennent les reactions hyperstatiques correspondant aux liaisons 
surabondantes rendent stationnaire l'energie interne. 



Cest le theoreme de Menabrea ou theoreme du travail minimum. 
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IV - 3 Exemples 



IV - 3.1 Poutre console 

Fleche en A et B? 



\f L/2 L/2 



B 



Rotation en A ? 



a) La poutre est sollicitee en flexion. Le moment flechissant est M 3 =-Px 

D'apres (12) 

r L Ml . r L P 2 x 2 . P 2 L 3 



w = r J^ dx = r 

Jo 9 FT Jo 



2EI 3 j0 2EI 3 



■dx 



6EL 



Done 



3W PL 3 



dP 3 EL 



(vers le bas) 



b) Pour determiner le deplacement du point B, on applique en B une charge fictive Q. 
PourO<x<L/2 M 3 =-Px et pour L/2<x<L M 3 =-Px-Q(x-L/2) 

d'ou 

W= (8P 2 + Q 2 +5PQ) 

48EL V J 



et finalement 
3W 



3Q 



L 3 



0=0 48EI 3 



(2Q + 5P) 



5PL 3 

Q=0 _ 48EI 3 



(vers le bas) 



c) Pour determiner la rotation de la section A, on applique en A un couple fictif M A . 
Pour 0<x<L M 3 =-Px-M A 

d'ou 



W 

2EI 3 

et finalement 

3W 



= fr^r(MA+Px) 2 dx = -i 



f p 2 L 3 



2EI 



3 V 



+ PM A L 2 +M A L 



0A = aM, 



PL 2 



u n 2EI 3 

M . =0 i 
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IV - 3.2 Treilli de barres 
E D 



->-lkN 



4 m 







3 4 m 




4 m 



On considere un treilli de barres articulees. 
Chaque barre est constitute du meme 
materiau et a la meme section. 
On cherche le deplacement du point C. 



2kN 



Q 



Pour ce faire, on applique en C deux charges fictives P et Q (avec P=Q=0) 
Equilibre du point C : cos45°F 3 =Q 



sin 45° F 3 + F 2 = P 



Equilibre du point D : F 4 = cos 45° F 3 + 1 



Equilibre du point B 



sin45°F 3 + F 5 = 
F 2 = cos45°F 6 +F 1 
sin 45° F 6 +F 5 = 2 



D'ou 



F, = P - 2Q - 2 , 



P-Q, F 3 =qV2, F 4 =l + Q, F=-Q, F 6 = V2(Q + 2) 



Par application de (10), 

W = Y F ' Li = — Yf 2 L. 



d'ou 



i=l 2S i E i 2SE i= , 



W = 2SE^ 2 +F > + ^ + F * +F * + ^ F ^ 



* Determination du deplacement vertical du point C 



3Q 



Soit 



p=o 2SE 

Q=0 



L 



f-4(P - 2Q - 2) - 2(P - Q) + + 2(1 + Q) + 4Q + 4^2(2 + Q) ] 

V Ap=< 





Q=0 



Ayc = ~si (5+4 ^ ) 



*Determination du deplacement horizontal du point C 



aw 



Ax r = _ 
c dP 



p=0 2SE 

q=o 



(2(P-2Q-2) + 2(P-Q)) 



p=o 

Q = 



Soit 



a 2L 
Ax r = 

SE 
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IV - 3.3 Application aux systemes hyperstatiques 



(0 



B 



La reaction en A est une reaction surabondante. 
q(x) = R A (x>"'-(o(x) 

d'ou 

x 2 

M, = R 4 x -a> — 



puis d'apres (12) 



W = f 

Jt 



l M 



2EI 3 



■dx 



2EI 3 V 



2 L n L coL ^ 
R; (0R A — + 

A 3 A 4 20 J 



Or comme le deplacement vertical du point A est nul 

' T3 



aw 

3R 4 



1 



d'ou 



R A =- 



2EI 3 V 



3wL 



L 4 ^ 

2R * — w = 



V - Equation de Bertrand de Fontviolant 



V - 1 Enonce 



L'equation de Bertrand de Fonviolant est une application directe du principe des puissances 
virtuelles. Nous avons vu (14), en considerant deux etats d'equilibre d'un meme systeme, que: 
f f f o 1 :!" da = f f f f 1 • u 11 da + f f F 1 ■ Q n ddQ. 

JJJQ JJJQ JJXl 

Si le systeme n'est soumis qu'a J forces ou couples ponctuels 
JJJ Q o:I n da = Jf/ ■ u\x\) + X C| ■ w u (x; ) 

i=l i=l 

d'oit, en developpant l'energie interne, on obtient: 



L |T,'T," f , 2 . ; | 





ES GS 



GS 



GI, 



EL 



+ ■ 



EL 



dX, = 



Xf/-u u (x;) + Xc:-6)"(x;) 



i=i 



i=l 



(16) 
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V - 2 Application: Evaluation des reactions hyperstatiques 
surabondantes 



CD 



^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ 



lAC L 




Iab 



B 



Reaction en B: Rb x , RBy 

Reaction en C: Mc,Rcx> Rcy 

Le systeme est hyperstatique de degre 2. 

r 



|R Bx +Rcx = o 



^R By +R Cy -03L=0 



G)L Z 



03 Rcy 



C Rex 



Rfiy 

B Rbx 



^M c -R By L + R Bx h+— = 
* On defini le systeme I par le portique isostatique associe (on supprime l'articulation en 



B). 



0) 



1 c 

i 



surAB M' = R Bx y 

surBC M 3 I =R Bx h-R By x+^- 



B 



Rbx 



* On defini le systeme II par : 



B 



surAB M3=F t y 

surBC M 3 U = F,h - F 2 x 

Le point B ne se deplacant pas 



En negligeant l'influence des effort tranchant et normal, par application de (16), les points 
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B et C etant fixes, on obtient 

3 3 dy + 
Jo pT J Jo 



a_3_j,. ■ f 3 3 dx = 



EI 



AC 



J — [(RB.yXF.y^y + 1 — , [R Bx h - R By x + — j fch - F 2 x) ,dx = 

A R .AC L I 



cox 2 \ 



R Rx F,h 2 



21 



+ • 



L 



L 



L 3 



L 4 ^ 



AB 



AC 



F,R Bx h z L - F]hR By — + F,h©- - F 2 R Bx h — + F 2 R By y - F 2 ®jj = 



I A cR Rx h 2 
2I AD 



L 2 



L 



+ R Bx h L - hR H „ — + h(0— | + F 



L 



L 3 



L 4 ^ 



R Bx h— + R R „— -(0— ; = 



^By 



Ceci doit etre vrai en particulier pour Fi=0 ou F2=0, d'ou 
I ^AcRjfcjL + 6R Bx h 2 L - 3hR n „L 2 + h(0L 3 = 



v By J 



■12R Bx hL 2 + 8R By L 3 - 3eoL 4 = 



_ . , I Ar h 

Soit, en posant k = ' 



R Bx - 



coL 



4h(4k + 3) 



et R R „ = 



3(0L(k + 1) 



By 2(4k + 3) 



V-3 



Application: Determination des emplacements et rotations 



Si en un point d'abscisse curviligne X l , on applique une force ponctuelle unitaire pour le 
systeme I, on obtient: 

J o ([-)dX 1 =u n (X 1 p ) 

ou u u (Xf) est le deplacement du point dans le sens d'application de la force unitaire 
Exemple: 

On cherche, dans le cas d'une poutre console chargee uniformement, les deplacements et 
rotations aux point A et B. 

CO 



L/2 



B 



L/2 



Pour le systeme courant, que nous appellerons le systeme II, on a: 
M" = -eo— 



* Fleche en A 

Soit le systeme I suivant: 
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Dans ce cas M3 = -1 x x . D'ou: 

^M^ dX lxu n 
Jo EL 1 



soit 



n 1 f L , x 3 COL 4 

u, = co — dx = 

A EI Jo 2 8 EL 



B 



* Fleche en B 

Soit le systeme I suivant: 



B 



Dans ce cas M l 3 = -1 x(x . D'ou: 

u!=- -f co— - (x — "2/ 1 dx = — — f co— (x-^)dx = 

PT Jo 9 z/ pT Jl/2 9 v 



17coL 4 
384EL 



* Rotation en A 

Soit le systeme I suivant: 



r- 



Dans ce cas M l = -1 . D 



oU: 



B 



0" = [ 3 3 dX, = [ co— dx = 

A Jo EL 1 EI J 2 



coL 3 
6EL 



m 



* Rotation en B 

Soit le systeme I suivant: 



¥■ 



B 



Dans ce cas M 3 r = -1 x(x D'ou: 

„ f LMiM? , 1 f L X 2 , 

9" = f 3 3 dX, = [ co— dx = 

B Jo pj 1 pi Jl/2 9 



EL 



7coL 3 
48EL 
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Flambement 



FLAMBEMENT 



I - Stabilite d'une poutre en compression 



On etudie le comportement d'une poutre droite soumise a une compression P et une charge 
repartie transversale q(x). On se place dans la configuration deformee de la poutre. 




D'apres (11-16) nous savons que l'equi 
Pourchaqueint ervallejxj , X J +1 



ibre se traduit par: 
f f +f , =0 



PourXji = l,...,J 



(a) 

c f +M 4 +ej aT = (b) 

F i +f(X 1 1 + )-f(X 1 1 ") = (c) 

C i +M(X 1 1 + )-M(X 1 1 ") = (d) 

Or dans le cas present, pour un tron§on, nous obtenons: 
f f =q(x)e 2 ,c f =6,C° =C 1 =6,F° =-F 1 =P 

et en considerant la deformee 
1 



U2,l 





SOlt 



Ti,i = 

T 2,i = -q( x ) 

M^i+^-TjU^ 
T 1 (0) = -P 







ou encore 



T { =-P 

T 2,i = -q( x ) 

M 3!ll -q(x) + Pu f 2ai =0 
et enfin en utilisant l'expression du moment flechissant (III -20) 



Mo = Elou 



f 



(1) 



(2) 



(3) 



(4) 



3 u 2,11 
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on peut conclure que 





EI 3 u 2,iiii+ Pu 2,n = q00 


(5) 


Si la charge laterale est nulle, l'equation d'equilibre se traduit par: 






3 dx 4 dx 2 


(6) 


Une solution triviale de (6) est u 2 = , mais dans le cas general on trouve: 




u 2 (x) = Cj +c 2 x + c 3 sin(nx) + c 4 cos(nx) 


(7) 


ou 


-I' 
V EI 3 


(8) 



Pour determiner les 4 constantes, nous utiliserons les conditions aux limites suivantes: 

* Le deplacement lateral de la fibre moyenne 

u 2 (x) = Cj +c 2 x + c 3 sin(nx) + c 4 cos(nx) (9a) 

* La pente de la deformee de la fibre moyenne 

u 21 (x) = c 2 +nc 3 cos(nx) -nc 4 sin(nx) (9b) 

* Le moment flechissant 

u 211 (x) = ^ 3 ^ X ~* = -n 2 c 3 sin(nx)-n 2 c 4 cos(nx) (9c) 
EI 3 

* L'effort tranchant (en utilisant (3)) 

T 2 (x) M 31 (x) Pu f 24 
EI 3 EI 3 EI 3 

= -u f 2411 -n 2 u 2a (9d) 
= n 3 c 3 cos(nx) -n 3 c 4 sin(nx) -n 2 (c 2 +nc 3 cos(nx) -nc 4 sin(nx)) 
= -n 2 c 2 



II - Etude de quelques cas simples 

II-l Colonne Rotule-Rotule 

Nous considerons une poutre rectiligne ayant une rotule a chacune de ses extremites. Nous 
recherchons alors la charge critique de flambement que peut supporter cette poutre. 
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En chaque rotule le emplacement lateral et le moment sont nuls. Done, si on elimine la 
solution triviale u 2 =0, l'equilibre se traduit par une deformee de type (7) avec pour 
conditions aux limites: 

u|(0) =Cj +c 4 =0 
M 3 (0)/EI 3 =-n 2 c 4 =0 
u 2 (L) = Cj +c 2 L + c 3 sin(nL) + c 4 cos(nL) = 
M 3 (L)/EI 3 = -n 2 c 3 sin(nL) -n 2 c 4 cos(nL) = 



soit, 



cj =0 
c 2 =0 

< 

c 3 sin(nL) = 
c 4 =0 

si on elimine la solution triviale, on obtient un equilibre pour nL=N7L C'est-a-dire, en 



utilisant (8), pour chaque compression telle que: 



P = 



(N7t) 2 EI 3 
L 2 



et la deformee est telle que: 
u 2 (x) = c 3 sin 



(10) 



(11) 



On remarquera que le deplacement lateral est indetermine! 
En conclusion: 

7l 2 EI 3 

* Pour P < P cr = l'equilibre est stable la poutre ne flambe pas. O 



■O 



7l 2 EI 

* Pour P = P cr = r-^- la poutre flambe 

L 2 

47l 2 EI 3 

* Pour P cr < P < — 2 - l'equilibre est instable 




47TEI, 

* Pour P = — — la poutre flambe 

L 2 




etc 
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II-2 Colonne Encastree-Libre 

Nous considerons une poutre rectiligne encastree a l'une de ses extremites. Nous 
recherchons alors la charge critique de flambement que peut supporter cette poutre. 




A l'encastrement, le deplacement lateral est nul, ainsi que la pente et l'effort tranchant. A 
l'extremite, la poutre est libre de tout moment. Done, si on elimine la solution triviale u 2 = , 
l'equilibre se traduit par une deformee de type (7) avec pour conditions aux limites: 
U2(0) = Cj +c 4 =0 
u 2,l(°) = c 2+ nc 3 =0 

T, /EL, =-n 2 c 9 =0 



soit, 



M 3 (L)/EI 3 = -n 2 c 3 sin(nL)-n 2 c 4 cos(nL) = 



c 2 =0 
c 3 =0 

c 4 cos(nL) = 



71 



si on elimine la solution triviale, on obtient un equilibre pour nL = (2N - 1) — 
La charge critique est done atteinte pour: 



Pcr = 



7t 2 EL 



et la deformee est telle que: 



u 2 (x) = c 1 



1-cos 



2L 



(12) 



(13) 



J J 



II-3 Colonne Encastree-Rotule 



Nous considerons une poutre rectiligne encastree a l'une de ses extremites, dont l'autre 
extremite est astreinte a reste dans l'axe. Nous recherchons alors la charge critique de 
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flambement que peut supporter cette poutre. 




A l'encastrement, le emplacement lateral et la pente sont nuls. A l'extremite, la poutre est 
libre de tout moment et ne peut se deplacer lateralement. Done, si on elimine la solution 
triviale u 2 = 0, l'equilibre se traduit par une deformee de type (7) avec pour conditions aux 



limites: 



u|(0) = Cj +c 4 =0 

u 2,i(0) = c 2 +nc 3 = 

u 2 (L) = Cj +c 2 L + c 3 sin(nL) + c 4 cos(nL) 

M 3 (L)/EI 3 = -n c 3 sin(nL)-n c 4 cos(nL) = 



soit, 



Cj = c 3 tg(nL) 

< c 2 =- c 3 n 

c 3 (tg(nL)-nL) = 
c 4 = -c 3 tg(nL) 

si on elimine la solution triviale, on obtient un equilibre pour nL= 1,4371. 
La charge critique est done atteinte pour: 



Pcr = 



7U 2 EI 3 



(0,7L) Z 

et la deformee est telle que 

f 

u 2 (x) = c 4 COS 



■0,223 sin 



1,43tcx^| x 



+ ■ 



(14) 



(15) 



III - Generalisation: Formule d'Euler 

Historiquement, on appelle "formule d'Euler" la charge critique d'une poutre Rotule- 
Rotule. 



Pour l'appliquer aux differentes combinaisons possibles d'appuis, on definit la charge 
critique comme etant: 
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Pour determiner le facteur K on peut effectuer les developpements mathematiques, ou bien 
utiliser des arguments de symetrie. 

Par exemple, dans le schema ci-apres on represente en noir la poutre etudiee, et en gris la 
symetrie utilisee pour retrouver la colonne rotule-rotule. 




K=2 





K=l 
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Autre ecriture: 



Pcr = 



7t 2 EI 3 
(KL) 2 



7C 2 ES 



(16) 



ou S est la section et r 
KL 



est le rayon de giration de la section. 



est appele le coefficient d'elancement de la colonne. 
r 

La contrainte critique correspondante est alors 



cr 



7l 2 E 

m 



(17) 



IV - Exemple 



On cherche la charge maximale que peut supporter la structure ci-apres. Ce systeme est 
constitue de deux colonnes (poutre en I) encastree a leur base, supportant une poutre supposee 
infiniment rigide. 




Comme l'indiquent les courbes pointillees, on peut identifier deux mecanismes de 
flambement, l'un dans la direction laterale et l'autre dans la direction arriere. 
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* Flambement vers l'arriere 

On se trouve dans la configuration d'une poutre encastree-libre, done K=2 et le coefficient 
d'elancement de la colonne dans la direction arriere est 



KL KL 2x4 10 3 




= 143,6 



* Flambement lateral 

On se trouve dans la configuration d'une poutre encastree, done K=l et le coefficient 
d'elancement de la colonne dans la direction laterale est 
KL KL 



r y 




= 98,5 



* Charge maximale 

Le coefficient d'elancement ayant une valeur maximale pour le flambement vers l'avant, il 
s'agit du cas le plus defavorable. On peut alors determiner la charge critique supportee par une 
colonne: 

7i 2 ES % 2 x 200 10 9 x 2960 10" 6 



P rr = - = — = 283343 N 

KLf 143,6 2 



V r J 

Comme chaque colonne supporte la moitie de la charge totale appliquee, on peut conclure 
que 

Pmax= 566,5 kN 
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COMPORTEMENT AU DELA DU 
DOMAINE ELASTIQUE 
CHARGES LIMITES 

(Chapitre en cours de redaction .... Approche de l'analyse limite par l'exemple) 

I - Introduction 

I-l Criteres de defaillance 

- Criteres cinematiques : Definis par le cahier des charges 

- Criteres de resistance : 

*Rupture statique brusque (comportement elastique fragile 

* Exces de deformation plastique 

* Rupture par fatigue 

* Fluage 

* Ecroulement 

* ... etc ... 



1-2 Comportement du materiau 




Comportement reel Comportement elastique-plastique parfait 
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ll - Analyse limite en traction 



II- 1 Analyse elastique 




Ce systeme est hyperstatique de degre 1 . 

En appliquant le Principe Fondamental de la 
Statique pour traduire l'equilibre du point A, on 
obtient : 



4i 



+ T AC + 



4i 

Tab j Tad _ q 



— F = 



V2 V2 




Pour resoudre ce systeme hyperstatique, on ecrit la compa- 
tibilite des deplacements en A. 



; < 



AL 



AC 



4i 



avec 



2 Lab ^ = T AB LV2 
ab ES 

soit en definitive 

T - F 
x ab - 



et AL AC = 



TacL 
ES 



T AC = 



2 + V2 
2F 

2 + V2 



La contrainte normale dans chaque barre est constante (etat de traction pure). Le premier 
barreau qui atteindra la limite d'elasticite est le barreau AC. On peut alors determiner la 
charge limite d'elasticite F e , ou encore la charge pour laquelle apparait de la plasticite: 



L AC 



soit 



F e = So e 



1 + 



J_ 



et le deplacement correspondant 
AL ACp 



o e L 



(2) 



(3) 
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II-2 Analyse elastique-plastique 




soit, 



Pour une charge F>Fe le barreau AC est 
entierement plastifie et ne peut supporter de 
charge supplementaire. L'effort normal dans le 
barreau AC est done constant et egal a Sa e . 

En appliquant le Principe Fondamental de la 
Statique pour traduire l'equilibre du point A 
(isostatique), on obtient : 

F = 



TAB +So e + TAD 



4i 



AB 



+ 



AD 



V2 V2 



= o 



(4) 



Lorsque le barreau AB (ou AD) atteint la limite d'elasticite, alors les trois barreaux sont 
plastifies et la structure n'est plus en etat de deformation limitee (mecanisme de ruine). 
On peut determiner la force limite pour laquelle on atteint ce mecanisme de ruine: 



On cherche la force Fi telle que AB 

S 



= G e . Soit: 



Fi =So e (l + V2) 



et le deplacement correspondant 

r 2 °e L 
AL Aq = ^ 2AL AB = ~ 



(5) 



(6) 
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Attention: Le emplacement n'est plus proportionnel a la charge appliquee. Le principe de 
superposition des contraintes ne peut plus etre utilise. 

II-3 Decharge 

On se place juste avant la charge ultime. Les barreaux AB et AD sont elastiques. Le 
barreau AC, par contre, est plastifie et a subi un allongement permanent. A la decharge, le 
comportement des trois barreaux est elastique. La barre AC sera comprimee et les barres AB 
et AD tendues. 

Pour trouver les efforts normaux residuels, on remplace F par -Fi dans (1) pour obtenir le 
retour elastique, puis on superpose avec (4) 
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III - Analyse limite a la torsion 



III-l Generalites 




On a vu que pour un cylindre soumis a de la torsion pure 
Mi 

t = Gi%i = — !-r 



Casl: On augmente le couple de torsion de telle sorte que les 
fibres exterieures atteignent la limite d'elasticite a e . 



Comme I, = 



nous obtenons done le moment limite 



elastique 



Mie =— — O e 



et 



e e =f 



M, 



-dXj = 



(8) 



(9) 




GIj 1 GR 

Cas2: On augmente le couple de torsion. La deformation 
augmente, done Tangle de torsion augmente, mais la contrainte 
ne peut depasser la valeur limite G e . 



Cas3: La limite d'elasticite est atteinte dans toute la section du 
cylindre. Toutes les fibres supportent done une contrainte a e et 

rien ne s'oppose a la rotation de la section. Le couple limite est 
alors: 

M 1L =jj(X 2 a 13 -X 3 a 21 )dX 2 dX3 

= jj (r cos 0a e cos0-rsin0(-sin9))rd0dr 

rR„2, 



27CO e j K r / dr 



M 1L = 



2 a 4 

3 e 3 le 



(10) 



II existe done une reserve de 33% avant l'ecoulement de la section ne soit complet. 
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III-2 Exemple 




Soit un cylindre plein de longueur L 
encastre aux deux extremites. on lui 
applique, a un tiers de sa longueur, un 
couple de torsion Mb. 



et 



Analyse elastique 

M A +M C =M B 

^ac = ^ab + ®bc = 



Soit <^ 



M A =-M B 
M c =^M B 



A la limite du comportement elastique, la section A commence a plastifier. Soit MA=Mi e 

et 

3 A/f 3ttR 3 
M Be =-Mi e =^— a e 

* Analyse limite 

Le seul mecanisme de mine possible est lorsque les deux sections en A et C on atteint la 
valeur limite. Soit: 

M A =M C =M 1L =^M le 

et 

M BL =| M le 

M BL /M Be =l,78 

Done, a partir de l'apparition de la plasticite, on dispose encore d'une reserve de 78% avant 
la ruine complete de la structure. 
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IV - Analyse limite a la flexion 



IV- 1 Generalites 




Dans le cas de la flexion plane, la contrainte varie 
lineairement avec l'epaisseur, 



o n =- 



M, 



I 



ainsi, bien sur, que la deformation axiale 

e ll = -~ = ~%3 X 2 

E 

Le moment maximum est atteint lorsque les fibres extremales subissent une contrainte Oe. 
Soit: 



M 3e =±G e 



I- 



(11) 



w=l3/h est appele le module de flexion. 

Si on depasse le moment M3e alors les fibres inferieures et superieures plastifient. Lorsque 
toute la section est plastifiee, on dit que Ton a une rotule plastique. 




Ge 

Dans ce cas: 

M 3L =-||X 2 o 11 dX 2 dX 3 =-o e {{X 2 dX 2 dX 3 

s s 

soit 

M 3L =2o e S x (12) 

ou Sx est le moment statique de la moitie de la section droite par rapport a l'axe X3. 

Plus generalement, on peut ecrire 
M 3L =KM 3e 

ou K est le facteur de forme de la section (dependant uniquement de la geometrie). 
Section rectangulaire pleine K=l,5 
Section circulaire pleine K=l,7 
Section en I K= 1,1 a 1,2 
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IV-2 Exemple: Methode "pas a pas" 



ZT 

A 



D 



LA 
B 



On etudie une poutre sur trois appuis 
supportant une charge F et on cherche le 
C mecanisme de ruine et la charge limite. 



* Etude elastique 




R 



F + R R +R r = 



- — F + LR B +2LR C =0 



q(x) = R A (x) 1 -F(x-L/2> 1 +R b (x-L) 1 
M 3 (x) = R A (x) 1 - F(x - L/ 2) 1 + R B (x - L) 1 

f 1 
EI 3 u 21 (x) = - 



R A (x) 2 -F(x-L/2) 2 +R b (x-L) 2 +Cj 



EI 3 u^(x) = ^[r a (x) 3 -F(x-L/2> 3 +R b (x-L) 3 +3C 1 x + C 2 



Pour determiner les constantes et l'inconnue hyperstatique, on utilise les conditions aux 
limites: 

u f 2 (0) = u f 2 (L) = u f 2 (2L) = 
et on obtient 



Ra = 



13F _ 22F 

, R B - 



Rr — , C 1 — , C 2 — 

32 32 2 



32 " 32 

En observant le graphe du moment flechissant ci-apres, on constate que le maximum est 

atteint dans la section D. 

M3. 13FL 




- 3FL 

32 
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Charges limites 



* Premier Pas: 

On augmente la charge F jusqu'a ce que la section D plastifie entierement. Dans ce cas on 



M l=— FiL 



13 
64 

soit la charge 

64 M L 



Fi = 



13 L 



Pour cette charge le moment flee his sant dans la section B est : 

Mo(L) = - — Mi 
3V 13 



* Deuxieme Pas: 



2 F 




Rotule 
plastique 



R'B 



Rc 



On augmente alors la charge Fi de AF. 
En decomposant le probleme, on se 
trouve maintenant dans le cas d'une etude 
elastique avec une rotule dans la section 
D. 



Les equations d'equilibre nous donnent les reactions isostatiques; 
R' B =-|AF, R' C =|AF 

et le diagramme des moments flechissants 
M3. 




- 2 FL 



Done, par superposition avec le probleme precedent, on constate que la section B plastifie 

ensuite lorsque M3(L)=-Ml, e'est-a-dire: 

f> 1 
M T AFL = -Mi 

13 2 

ou encore pour un accroissement de charge 
13 L 

Lorsque les section D et B sont plastifiees, la structure s'effondre. On a atteint alors la 
charge limite 



F L = Fj + AF = 6 



M, 
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La fleche dans la section D est 

u f (L/2) = ^L^ + ^ = AMLtl 

1536 EI 3 8EI3 24 EI 3 



IV-3 Theoreme energetique 

D'apres le principe des travaux virtuels, la somme du travail des efforts exterieurs et de 
l'energie interne est nulle. Si on se place dans le cas d'une structure soumise en flexion a un 
etat limite (mecanisme de mine) alors seules les rotules plastiques contribuent a l'energie 
elastique. 

Soit une structure soumise a N forces d'intensite Pi (l<i<N) et presentant un mecanisme de 
ruine a M rotules plastiques. On note 0j (l<j<M) la rotation d'une rotule plastique et 8i le 
deplacement au point d'application d'une force. 

On peut alors ecrire: 

N M 

i=l j=l 



Exemple d'application: 

CO 

^^^^^^^^^^^^ 



A 



B 



C 




Le mecanisme de ruine apparait lorsque les trois 
sections A,B et C plastifient entierement. Dans ce cas la 
deformee est telle que: 

Travail des forces exterieures 

L L/2 L/2 1} 

W e = }cou2(x)dx = 2 j coxtgGdx = 2 {coxGdx = Geo — 
4 

Travail des forces interieures 
-Wj =M A G A +M B G B +M C G C =-M A G + 2M B G + M c G 



Or a 

M A =-M L , M B =M L , M c = M L , done 



l'etat 



limite 



on 



Geo 4M T G = 

4 L 

et la charge limite est done: 
M T 
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Charges limites 



Verification par la methode pas a pas: 
Ra 

A 

CO 




Le systeme est hyperstatique. Par symetrie on 
a naturellement: 

R -R - &L 

K A - K B - ~7T~ 



M A =-M B 

Afin de determiner l'inconnue hyperstatique nous cherchons la deformee: 
q(x) = -M A )x(- 2 +R A )x(- 1 -C0>x(° 

M 3 (x) = -M A )x(° + R A )x( 1 -|)x( 2 



EI 3 U2!(x) = -M A x + - 



R 



EI 3 u t 2 (x) = -^ 



M 



■x +■ 



2 
R 



2 CO 3 

-x x 3 

6 

A x 3 --^x 4 
24 



2 6 

La condition de symetrie u 21 (L/2)=0 nous amene a 



M A = 



12 



et 



M 3 (x) = -| 




En considerant le comportement elastique 
jusqu'a la rotule plastique, les premieres 
sections a plastifier sont les sections 
extremites A et C, pour une charge COe telle 
que : 

|M 3 (0)| = |M 3 (L)| = M L 
soit 



co e = 



12M, 



La fleche maximale est 



co e L^ 



M L L Z 



384EI 3 32EI 3 

Dans une deuxieme etape, on se place dans le cas ou la poutre est soumise a une charge 
repartie, et du fait des rotules plastiques aux extremites, soumise a 2 couples sur les sections 
extremes en appui simple. 
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co 




Ce systeme est isostatique. Le moment 
^ ^-^Ml flechissant est: 

AT / N A/T ^ 

M, x =-M, + x 

3W L 2 2 

Le maximum est naturellement atteint dans la 



section milieu B: 

M 3 (L/2) =-M L + 



G)L 2 (OL 2 odL 2 



8 8 



-M, 



On aura un mecanisme de ruine quand M3(L/2)=Ml, soit 
16M T 
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Annexe 



ANNEXE 
FORMULES ESSENTIELLES EN 
MECANIQUE DES MILIEUX CONTINUS 



1. Coordonnees cartesiennes orthonormees 



OM = xe x + ye y + ze z 

* Soit v = v x e x + v y e y + v 7 e 7 un vecteur, alors 



grad(v) = Vv = — -f-Bj ® e = = 



3v x 


8v x 




dx 


dy 


dz 


3v y 


dv y 




~dT 


~W 


~aT 


8v z 


dv z 


8v z 


dx 


dy 


dz 



et 



divv = = Tr(grad(v)) = ^ + + ^ 
dxj dx dy oz 



Av = div(grad(v)) = 



3 Vj 

axjaxj 



ej = Av x e x + Av y e y + Av z e 



* Soit f une fonction scalaire, alors 



grad(f) = Vf =m 
dx: 



K 

dx 

K 

dy 

K 

dz J 



* Soit T = T ij e i ®ej 



T T T 

A xx A xy A xz 

T T T 

y x yy y z 

T T T 

A zx A zy A zz 



et Af = div(grad(f )) = — — — = —- + — - + — - 

dXjdxj 3x 2 dy 1 dz l 

un tenseur symetrique du deuxieme ordre, alors: 



= dT„ 
div(T) = -^e 1= < 
3 Xj 



3T VV . 



■ + ■ 



xy 



dy 



dx 

9T yx r)'i; 



+ 



dT 



dx 



dz 

_l yy _|_ yz 



dy 



9T, V + ^zy_ + ax 



9x 



8y 



dz 



= 3 Tj: 

et AT = — ej ® e , = 

dx k dx k 



AT XX AT xy AT XZ 
AT yx AT yy AT yz 
AT ZX AT zy AT ZZ 
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2. Coordonnees cylindriques 



OM = re r + z5 7 et 

z 3r 



30M _ 130M „ 30M 



30 



e 6' 



3z 



d(OM) = 5 r dr + rd0e e + e z dz 



35 



3r 
35 



^ = 



35 



e 



36 

35 r 



o 



3r 
35q 

30 

35 A 



= 



-e. 







35_ 2 

3r 
35 



= 



30 

35, 



^ = 







3z 3z 3z 

* Soit v = v r 5 r + v e 5 9 + v z 5 z un vecteur, alors 



grad(v) = Vv 



3v r 
3r 



1 



r \dQ 
1 



— v e, 

1 dv z 

r ae 
3v r 



3^ 
9z 



et 



divv = Tr(grad(v)) = — ^ + ^ + - 

or r r o0 



1 3v fi 3v, 



+ ■ 



Av = div(grad(v)) = 



Av r - 



2 3v e ^ v r ^ 
r 2 30 r 2 y 



3z 



e r + 



Av e + 2 



2 3v, 



30 



ve 
r 2 J 



5 r +Av z 5 z 



* Soit f une fonction scalaire, alors 

grad(f) = f e r +I^e e + f 5 Z et Af = ^ + + J-^l + ^L 
9r r r 99 6 9z z 3r 2 r 3r r 2 o0 2 oz 2 



* Soit T = T ij 5 i ®5j 



div(T) = 



3r 



T 

zr 

1 9T r6 3T rz 



x zz 



un tenseur symetrique du deuxieme ordre, alors: 



+ - 



3T, 



3r 
9r 



+ 



r ae 

1 3T, 



+ 



3T, 



6z 



+ 



2X, 



9z 

z6 I d T zz I T zr 
99 dz r 



3. Coordonnees spheriques 

_ 30M _ 130M _ 1 30M 



OM = r5 r et^— = e r , — — = 5 e ,— -— — = 5 
3r r 30 rsin0 3cp ^ 

d(OM) = 5 r dr + rd05 e + e^r sin 0dcp 
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Annexe 



35 r 

~3r~ 

35 r 

~ae 

35 r 
3cp 





g e 
sin 05 



35 



e 



3r 
35 



= 



35 



3r 
35 



= 



<P _ 



<P 



30 
35 fl 

— = COS0e<p 



30 

35 



= 



3cp 



— = sin05 r -cos05 e 



* Soit v = v r 5 r + v e 5 e + VjpSjp un vecteur, alors 

3r r Vae 8/ r \sin9 3(p >/ 



grad(v) = 



r VsinO 3(p 
9v e 1 Pv fl ) 1 ( 1 3v e 
3r r V30 r VsinO 



3(p 



cot g6v 



dv 9 
3r 



1 9v <P 

r ae 



ll i 3v <p 



r I sin 9 3(p 



+ cot s9v fi + v. 



et 



divv = Tr(grad(v)) = ^ + 2^ + 1^- + - J-^L + cot g0^> 
3r r r 30 rsin0 3cp r 



Av = 



Av r - 



1 3(v e sin6) + 1 9v 9 
in 6 39 sin 9 3q> 



Av M 9v r v e cose av 9 

6 r 2 l 36 2sin 2 9 sin 2 6 3cp 



Av ( „+ 



<P ' 2 • 



2 fa Vr 



+cotg0 



3va 



r-sin9v3(p 3<P 2sin9 

* Soit f une fonction scalaire, alors 



grad(f ) = 



K 

3r 
131 

r ae 

1 3f 



rsin9 3(p 



et Af = 



3 2 f 2 3f 1 3 2 f 1 3f 1 3 2 f 



• + + - 



3r 2 r 3r r 2 30 2 r 2 



+ — cotg0— + ■ 



Soit T = T ij 5 i ®5j 



div(T) = 



d T rr | 1 9T r6 



+ 



Tit T r g 


T 


T e r T ee 


T 69 


T(pr T,pg 


T 


1 9 T r<p 


+ 1 

r 


r sin 6 39 



30 r 2 sin 2 3(p 2 
un tenseur symetrique du deuxieme ordre, alors: 



2 Trr-Tee-T w +Tre co U 



"3T + 7 IT + TshTe + 7 Ll T ee-T w jcot g e + 3T r6 J 

r)T,„,„ , f j 



3T, 



3T m 



* + 1 ^ + i ^ + j K cot g9 + 3T, 



3r r 36 rsin9 3(p 



rep. 



4. Comment retrouver les formules en coordonnees 
cylindriques 



On note V = v r 5 r + v e 5 e + v z 5 z = Vjej avec i=r, 0, z et ,i = 



3r ' r 30 ' 3z 



Done, avec cette convention 5 r e = — et e e e = — - 



Chercher le gradient d'un tenseur consiste a augmenter l'ordre de ce tenseur, soit 



I.S.I.T.V. 



Resistance des Materiaux 



-87- 



grad(**) = (**) j ®ej 

Si on applique cette remarque a un vecteur, on obtient: 
grad(V) = (v i e i ) J ®e j 

En n'oubliant pas de deriver les vecteurs de base car nous sommes dans un systeme de 
coordonnees cylindrique, 

grad(V) = (viBj )j ® e j = VyBj ® Bj + v^y ®e ] = VyBj ® Bj + v^e ® e e 
= Vi jSi ® ej + v r B r6 ® e e + v e B e e ® e e 

V Vo 

= v ; je ; ®ej +— e e ®e e -8 r ®e 9 

Pour obtenir l'operateur divergence, il suffit de prendre la trace du gradient, 
div(**) = Tr(grad(**)) 



3v r 1 3v e 3v 



+■ 



z + ^ 



soit dans le cas d'un vecteur: 

div(V) = Tr(grad(V))=v H +^ = 

r or r o0 oz r 

Appliquons maintenant cette methodologie a un tenseur d'ordre 2. 

grad(T) =(T 1J H i ®H J ) k ®e k 

= T ij;k e i ®5j ®e k +T^ k ®e } ®e k +T ij e 1 ®e j>k ®e k 

= T ij;k e i ®Hj ®e k +T y e i>8 ®5j ®e e + 1^ ®e jje ®B e 

_ T r j ^ _^ Tq: _^ 

= Ty^e; ®ej ®e k +— e e ®ej ®e e -e r ®ej ® e e 

+ ^ L e i ®e e ®e e -— Bj ®e r ®e e 
r r 



Pour obtenir la trace de ce tenseur d'ordre 3 on contracte les deux derniers indices: 
divfrU Trfgrad(T)) = T ;i , 8: + ^-B fl - ^-B r + — e, 



^3T. 



e 9 — e r +■ 

r r r 

| l 3T r6 | 3T rz T ee | T 

3r r 30 3z r r 



+ 



+ 



V ui i OZ 

^T 9r | 1 3T ee | 3T, 
3r r 30 



9z 



L r9 



A 3T 



+ — — + 
3z r r 



9r b 



e e 



+ - 



l3T z6 3T 



3r r 30 



+ 



3z 



T ^ 

lz _|_ ^zr 



On peut done maintenant retrouver l'operateur Laplacien d'un vecteur 
Av = div(gradv) 



= v i,jj e i — — e e 



Av r 



r 

2 3v 



v e,e- , 
e r + 



r,e 



v 



e e 



v e,e + 



v 



— -e r H ej 



e 



30 



r 2 y 



e r + 



A , 2 ov r v e 
v 6 r 2 30 r 2 j 



e e +(Av e )e 2 
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